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INTRODUCTION 

Ip  this  report  I  have  attempted  to  sort  through  the  workings  of  a 
Synthetic-Aperture  Radar  (SAR).  The  motivation  is  to  define  those  things  which  are 
based  on  physics  and  sorting  out  those  which  are  engineering,  which  may  have  come  to 
be  conside.ed  basic  to  the  SAR  for  whatever  reason.  The  reason  for  this  attempt  was 
my  ineptitude,  in  a  discussion  with  an  engineer  who  is  new  to  the  program  to  explain 
the  motivation  for  such  a  program.  The  subject  related  to  optical  back-projection 
processing,  which  may  find  use  in  the  processing  of  polar  formatted  SAR  data. 

However,  our  discussions  centered  more  on  polar  formatting  and  its  whys  than  on 
optical  back-projection  and  its  attributes.  This  should  be  the  real  thrust  of  any 
optical  back-projection  program  at  this  point.  Application  can  come  later! 

I  have  felt  for  some  time  that  descriptions  of  SAR  have  suffered  from  the  same 
syndrome,  i.e.,  topics  which  are  really  engineering  conveniences  have  come  to  be 
treated  as  basics,  and  as  such  may  impede  attempts  at  alternates  because  of  fear  of 
the  impact  they  may  have  on  the  SAR.  The  approach  here  is  to  define  the  SAR  in  terms 
o*  its  most  basic  physical  requirements  and  then  explore  the  attributes  of  various 
existing  implementations  of  these  basics.  Perhaps  at  some  time  in  the  future  a 
reader  of  this  report  may  ask,  "That  being  the  case,  and  if  I  don't  violate  the  basic 
concept,  what  if...?"  There  is  much  that  can  still  be  done. 

In  the  first  sections  of  the  report  the  basic  physical  principles  of  SAR  are 
detailed.  The  radar  is  nothing  more  than  an  interferometer,  with  multiwave  length 
capability,  which  can  be  used  to  measure  the  optical  path  lengths  (optical  path 
differences)  to  scatter  points  in  the  radar's  field  of  view.  Data  on  the  target 
field  is  collected  from  various  aspects.  The  composite  of  optical  path  da^a  may  then 
be  considered  as  a  whole  and  we  have  shown  that  is  is  ideally  processed  with  a 
matched-filter.  The  physical  attributes  of  the  ideally  matched-filter  provide  the 
basic  limit  of  the  radar's  performance.  The  remaining  sections  of  the  report  roughly 
follow  the  historical  developments  of  SAR,  and  as  we  shall  see,  are  based  on 
engineering  approaches  taken  to  implement  the  desired  matched-filter,  and  are  not 
basic  radar  principles.  One  section  relates  the  synthetic  array  to  its  real  array 
analog. 
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The  data  of  the  earliest  SARs  were  processed,  using  optical  techniques  after  the 
data  were  literally  arrayed  in  a  range-azimuth  format.  A  reader  may  more  typically 
see  this  referred  to  as  a  range-doppl er  format.  I  avoid  this  designation  as  a 
Doppler  shift  is  not  a  basic  requirement  of  SAR.  The  characteristics  of  an  optical 
processing  implementation  and  its  digital  counterpart  are  detailed,  and 
formatting-based  physical  limitations  are  described  in  a  section  entitled  Processing 
Considerations.  The  drive  for  improved  resolution  lead  to  the  intermixed 
(range-azimuth)  matched-filter  approach  exemplified  by  polar  formatting.  A  limiting 
form  of  the  data  formatted  in  this  fashion  is  spatial  frequency  monotones  which  can 
be  processed  using  Fourier-transform-like  operations.  Subsequent  sections  of  the 
report  center  around  techniques  for  the  digital  implementations  of  these 
Fourier-transform-1 ike  processes.  The  discrete  Fourier-transform  and  the  efficient 
FFT  implementation  are  explained.  Data  interpolations  necessary  to  allow  use  of  fho 
FFT  algorithm  are  outlined.  A  final  section  deals  with  the  general  characteristics 
of  back-projection  processing.  Its  characteristics  are  illustrated  by  tying  it  hack 
to  the  basic  discrete  Fourier-transform  form. 

The  appendix  iltustrates  digital  processes  which  may  be  used  to  break  the 
discrete  Fourier-transform  operation  into  smaller  steps.  Keep  in  mind  that  all  the 
processing  discussions  pertain  to  various  approaches  to  efficiently  implement  a 
matched  filter  and  are  not  basic  to  SAR.  However,  the  attributes  of  any  one 
implementation  may  impact  the  performance  of  the  radar  implemented  in  this  manner. 
Throughout  the  report,  an  attempt  is  made  to  estimate  the  digital  computational 
complexity  of  the  various  matched-filter  implementations. 

To  keep  maximum  emphasis  of  the  basics,  we  have  assumed  a  linear  system  so  that 
the  analysis  can  proceed  by  considering  only  a  single-point  target.  We  noted  a  need 
for  the  collection  platform  reference  system  and  motion  compensation,  and  they  are 
ideal.  Inphase  and  quadrative  channels  are  needed  to  sort  out  ambiguous  signals,  as 
well  as  accommodate  arbitrary  scatter  phases.  We  explored  the  characteristics  of 
only  one  channel  (both  channels  contain  data  with  the  same  form).  Finally,  wo 
ignored  all  the  overhead  functions  which  have  to  occur  between  the  collection  of  the 
raw  data  and  its  appearance  as  a  sample  to  a  matched-filter  emulafor. 
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MICROWAVE  INTERFEROMETRY 

In  this  section,  I  will  describe  the  hardware  aspects  of  a  Synthetic-Aperture 
Radar  which  I  feel  is  its  most  basic  physical  form.  All  other  descriptions  can  bo 
derived  from  this  basic  form.  This  form  is  the  radar  which  functions  as  an 
i nterferometer  with  multiple  frequency  capability.  The  radar  collects  data  in  the 
form  of  a  differential  phase  between  a  received  signal  and  a  locally  established  and 
controlled  reference.  We  will  return  to  the  reference  form  below.  We  will  borrow 
■from  optical  interferometry  and  refer  to  this  differential  phase  as  the  optical -path 
difference  since  it  involves  1 ine-of-sight  distances  to  the  target's  positions. 

If  we  assume  that  the  radar  operates  in  a  linear  region,  we  can  view  the  received 
signal  as  the  composite  of  returns  from  many  scatter  points  in  the  radar's  field  of 
view.  Thus,  our  interferometer  will  be  simultaneously  collecting  the  optical-path 
difference  data  for  a  large  number  of  scatter  positions.  With  this  interferometer  wp 
can  systematically  collect  the  optical-path  difference  data  for  a  band  of  frequencies 
and  from  different  viewing  aspects.  The  totality  of  data  will  contain  optical-path 
data  on  each  scatter  point,  but  the  history  or  variation  in  the  optical-path  differ¬ 
ence  with  frequency  and  aspect  will  be  unique  for  each  target  position.  This 
uniqueness  of  optical-path  difference  history  with  frequency  and  aspect  anale  allows 
a  single  history  to  be  sorted  out  from  the  jumble  of  signals  available,  i.e., 
detected.  The  precision  of  this  sorting  out  process  is  limited  by  the  bandwidth 
transmitted  and  the  change  in  viewing  aspect.  Following  sections  of  this  report  are 
devoted  to  various  implementations  which  can  be  used  in  this  sorting  out  process, 
starting  with  the  ideal  and  its  limits.  The  limits  of  the  ideal  processing  scheme 
establish  the  basic  limits  of  the  physical  process  employed. 

Before  going  on  to  the  processing,  go  back  to  our  basic  radar  and  observe  that  wp 
require  not  only  coherency  between  the  signal  and  reference  at  the  time  of  any 
measurement,  but  an  absol ute  coherence  for  the  duration  of  the  total  observation.  If 
the  basic  frequency  of  the  radar  is  uncertain,  an  uncertainty  in  the  optical -path 
difference  histories  and  the  sorting  out  capability  may  be  negated,  or  at  best  be  in 
error  or  ambiguous.  This  absolute  coherency  is  established  by  tying  the  transmitter 
to  the  local  oscillator,  and  the  phasings  of  any  modulations  on  each,  to  a  local 
"absolutely"  coherent  source. 
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The  reference  for  the  interferometer  is  the  radar's  local  oscillator,  in  our 
case,  and  its  phase  can  be  controlled  to  provide  data  with  a  structure  consistent 
with  a  chosen  processing  scheme.  We  will  return  to  some  of  these  processinq  forms 
below.  For  now  we  can  assume  the  local  oscillator  is  locked  to  the  coherent,  source 
with  the  consequence  that  the  interferometer  is  measuring  the  total  phase  in  the 
two-way  path  between  the  measurement  position  and  the  various  scatter  points.  With 
the  interferometer  measuring  the  two-way  optical  path,  it  becomes  necessary  to  have 
an  instantaneous  knowledoe  of  the  measurement  position  as  the  various  frequency  and 
angle  measurements  are  made.  If  this  is  not  the  case,  deviations  in  measurement 
position  data  (absolute  path-length  change)  would  be  superimposed  on  the  desired 
data,  and  separation  may  be  impossible.  The  ideal  of  no  rel ati ve-range  change  may  be 
measurements  at  constant  range  and  depression  angle  while  viewing  at  various  aspects, 
i.e.,  a  constant  altitude  on  a  conical  surface  with  its  apex  in  the  target  area.  The 
turntable  SAR  emulator  is  such  an  implementation.  In  practice,  it  is  rather 
impractical  to  move  on  such  a  conical  surface.  Thus,  we  provided  an  inertial 
reference  system  as  part  of  the  radar's  hardware  implementation.  This  reference 
system  is  used  to  provide  an  instantaneous  knowledge  of  the  measurement  platforms 
position.  Inertial ly  derived  position  data  may  then  be  used  to  alter  the  measured 
data  so  that  it  appears  to  have  been  measured  at  some  other  position,  i.e.,  const,  nt. 
ranges  on  a  conical  surface  with  its  apex  in  target  space.  It  should  be  obvious  that 
this  modification  (addition  or  subtraction  of  phase  to  the  measured  optical-path 
difference)  also  reouires  a  coherent  system. 

The  hardware  described  above  is  the  totality  of  those  necessary  to  implement-  a 
Syntheti c-Aperture  Radar.  A  coherent  standard  is  required  as  is  a  pos i t ion-sens  inn 
system  if  an  absolute  measurement  position  cannot  be  assured.  The  standard  is 
ideally  coherent  over  a  time  long,  relative  to  the  total  measurement  time.  The 
position-sensing  apparatus  accuracy  must  be  some  small  fraction  of  the  operating 
wavel enqth . 

There  is  no  requirement  for  Doppler,  nor  is  there  any  restriction  on  the 
modulation  format.  The  Doppler  description  probably  arose  because  of  the  aircraft 
motion  and  its  convenience  in  specification  of  a  PRF.  The  PRF  is  a  samplino 
requirement  reflecting  the  rate  of  change  of  the  optical  paths  with  anglp,  which  is 
directly  relat.able  to  aircraft,  velocity  if  this  velocity  is  the  source  of  the  r ni'ir 
change.  The  synthetic-aperture  data  may  be  collected  from  a  stationary  pi  .«♦  11 
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which  moves  in  any  manner  between  measurements.  The  waveform  is  immaterial  cn  lorn 
as  its  frequency  components  can  be  referenced  to  the  same  frequency  components  of  the 
local  oscillator.  To  eliminate  ambiguities,  there  is  a  requirement  on  samplinq 
density  in  both  angle  and  frequency. 
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MATCHED  FILTERING 

Synthetic-Aperture  Radar  is  based  on  matched-filtering  concepts  but  for  a 
different  reason  than  those  used  to  formulate  the  matched-filtering. 

The  optimization  of  signal-to-noise  ratio  occurs  in  SAR  as  in  any  use  of  matched- 
filtering,  but  its  application  in  SAR  is  tied  more  to  the  properties  of  the  signal 
produced  by  the  filter.  Any  signal-to-noise  enhancement  is,  at  best,  related  to  an 
individual  resolution  element.  The  signal-to-noise  ratio  may  be  used  to  estimate  the 
performance  of  a  radar  against  a  target  of  reference  cross  section  and  as  such  can  be 
used  to  estimate  the  required  transmitter  power.  However,  these  S/N  considerations 
are  secondary  and  a  matched -fi 1  ter  would  probably  be  used  to  process  SAR  data  even  in 
the  absence  of  noise.  The  rationale  for  these  statements  will  be  illustrated  in 
the  remaining  paragraphs  of  this  section.  We  can  begin  by  looking  at  the  properties 
of  a  matched-filter  and  its  performance  when  working  with  the  optical-path  difference 
data  generated  by  our  coherent  microwave  interferometer. 

The  matched-filter  concept  was  originally  formulated  to  provide  a  linear  filter 
for  the  maximization  of  the  detected  signal-to-noise  ratio  when  a  signal  is  buried  in 
noise,  the  noise  being  characterized  as  samples  from  a  wide-sense  stationary  random 
process  [I],  The  knowledge  that  a  signal  does  exist  removes  the  constraint  of 
preserving  the  signal  form.  Under  the  assumption  of  white  noise  fa  flat  power- 
density  spectrum),  an  optimum  filter  has  an  impulse  response  form 

h  (  T  )  =  1  S(t  -  t)  ,  0  <_  T  T 

N  1  ( 1  ) 

where  S(t)  is  the  input  signal  and  T  is  its  extent.  A  filter  with  this 
characteristic  is  labelled  a  "matched-filter."  It  is  obvious  that  the  structure  of 
S(t)  has  to  be  known  to  define  the  filter.  The  signal  produced  by  such  a  filter  is 
represented  by  the  convolution  integral 
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where  a  fi'rre  dplay  has  bppn  included,  i.e.,  a  propagation  path  de1av  in  *ho  r-< 
a  radar  or  communi  cat  ions  system. 

The  properties  of  a  matched-filter  can  bp  i 1 1  ust  rated  hv  r  i  r, q  a-  ’e*f  • 

•-  t  ,  i.e.,  we  seek  some  knowledge  on  the  signal  in  the  t  ’me  slot  wher  v.e  m,i ,  e«g  eel 

a  signal.  He  should  recognize  that,  in  principle,  we  require  a  separate  t’’ter  cnr 

each  time  slot  of  interest.  This  requirement  is  what  drives  us  to  lock  onto  a 

communications  signal  and  provide  huge  amounts  of  processing  for  rada*-  ►singe 

compression.  With  t.  =  t,  and  the  substitution  t '  =  t-t  ,  we  find  (?'•  taking  *he  t,-,. 

Id  d 

SQ(t '  )  ,  S(td-  )  5 ( t '  -  - ) d t 


The  above  is  the  autocorrelation  function  for  the  signal  referenced  tp  the  evp<.  ted 
t me  delay.  At  t'  =  t^  the  indicated  integration  provides  an  output  proper*  i < •»  -  ’  < 
xhe  energy  in  the  signal  (volts  squared  times  time  or  current  squared  f’mes  ft'-e'.  . 
Had  we  been  in  error  in  estimating  ,  the  correlation  time  would  shift  in  t’nr  hv 
estimation  error.  We  also  know  the  autocorrelation  function  of  the  signal  ;e  relate 
to  its  power  spectral  density  by  the  Fourier-transform  relationship  •'  T T 

i?-.  • 

P  •  : )  =  S(v)e  d 

The  variable  v  may  be  hertz  or  ir verse  distance.  Drawing  upon  transform  proiiert  ii 
and  considering  a  fairly  uniform  power  densi ty  spectrum,  note  that  R(  1  }  and,  o^  F 
( t 1 1  will  display  significant  values  only  durina  a  time  interval  of  extent  VB 
centered  on  t1  -  0,  i.e.,  about  t  =  t^.  This  time/space  interval  will  contair  -rst 
of  the  energy  contained  in  the  original  signal.  The  amount  of  energy  regaining  in 
near-by  time/space  positions  are  side  lobe  levels  and  of  concern  when  wide  dynamic 
ranges  arc  dictated.  However,  these  side  lobes  are  not  has'o  to  an  understanding  c.'r 
',AP  and  for  purposes  of  illustration  we  can  assume  that  our  f  i  <  *er  comprccc-c 1  a  1 1  R 
incider',  s tonal  energy  into  a  time/space  interval  of  extent  ]'R,  where  R  is  the 
bandwidth  of  the  incident  signal. 


forsidering  ar  input  signal  of  fairly  uniform  power  we  may  reure'ent  tfu 
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We  require  a  minimum  bandwidth  of  B  to  pass  the  input  signal.  In  the  correlation 
time  window  the  output  signal-to-noise  ratio  becomes 


(St)  (1/B) 
NB 


The  noise  power  is  uncorrelated  and  remains  uniformly  spread  over  the  interval  T 
whereas  the  signal  energy  is  compressed  into  a  single  time  window  of  extent  (1/B). 
This  gain  in  signal-to-noise  ratio  in  this  time  slot  comes  at  the  expense  of  the 
signal-to-noise  ratio  at  all  other  times  in  the  original  signal  extent  T.  While  thi'- 
is  of  no  consequence  for  detection,  it  is  unacceptable  if  the  original  signal 
structure  has  to  be  preserved.  For  the  purposes  of  SAR  we  only  need  to  note  that  a 
properly  designed  matched-filter  compresses  a  signal  of  extent  T  into  one  of  extent 
1/B  ,  i.e.,  a  compression  equal  to  the  time/space-bandwidth  product  of  the  signal. 

The  time-bandwidth  product  of  a  monochromatic  pulse  of  length  T  approximates  1  and  no 
compression  is  possible. 


Let  us  now  adapt  the  above  to  this  microwave  interferometer  we  call  a  SAR.  We 
can  invoke  our  linear  system  arguments  and  consider  the  signal  generated  by  a  sinole 
frequency  at  a  single  aspect.  This  signal  is,  in  general,  complex  and  may  be 
represented  with  the  phasor  form 


where  is  the  optical  path  length  including  array  scattering  phase  and 

is  the  reference  phase.  To  avoid  tracking  in-phase  and  quadrative  channels  let 
us  assume  a  scattering  phase  of  zero  and  work  with  only  the  in-phase  channel.  We 
note,  however,  that  in  any  implementation  we  will  require  both  in-phase  and 
quadrative  channels  with  identical  forms.  The  matched- fil ter  for  a  signal  with  *he 


^  ipj  yj  v.  y;  if  1  y  j  y ;y  ■  ■  ■  ■  » 
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fomi  of  (7)  is  approximated  by  the  complex  conjugate  of  the  signal.  Our  proressing 
is  the  detection  of  a  signal  from  a  particular  point  in  radar  space,  with  the  sicca; 
buried  in  a  composite  of  signals  of  the  same  complex  form.  We  merely  weiqht  thr 
composite  of  each  measurement  with  the  complex  conjugate  of  (7),  record  the  value, 
and  then  sum  the  results  achieved  when  the  same  weighting  is  applied  to  all 
measurements.  Each  scatter  point  will  have  its  unique  matched  filter. 


To  illustrate  the  azimuth  resolution  of  a  SAR  we  can  consider  we  are  gathering 
data  on  the  returns  from  a  point  P,  as  a  single  frequency  measurement  system  tains 
various  positions  along  a  straight  line  such  as  shown  in  Figure  1.  In  the  case  of  an 
actual  radar  we  will  motion  compensate  the  data  to  make  it  appear  as  though  it  had 
been  collected  while  the  aircraft  was  flying  along  a  straight  line.  Our 
interferometer  is  measuring  the  optical  path  v  with  the  signal  of  interest 
approximating 


S 


Ae-2kr 


t  '• 

'  ,  ->  J 


2  ?  u 

here  v  =  (R  +  x  )  w  is  the  radar  radian  frequency  2Hv  ,  and  c  is  the 

velocity  of  light  in  medium  of  interest.  The  spatial  frequency  of  this  signal  ai  ary 

position  along  the  collection  aperture  is 

1  „  v  x _ 

vs  2tt  3v  ~  "  c  ^2  +  R2^1/2 

2  v 

=  -  —  sin  radians /me ter 

For  small  9  we  see  that  the  spatial  frequencies  are  a  fairly  uniform  function  of 
"*  i.e.,  have  a  fairly  uniform  power  density  spectrum.  In  an  interval  a e  as  shown  ui 
Figure  1,  we  will  observe  spatial  frequencies  in  the  range 


sin  o  radiaAS  ,  v  <  2v  in  radja_ns 
H  m  s  —  c  H  m 


'  1  "J  'i 


With  a  matched  filter  we  can  then  anticipate  an  isolation  of  the  sianal  sour*  e  to  < 
position  uncertainty  (azimuth  resolution) 


9 


B 


Collection  Path 


tic -Aperture  Data  Collection 
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Ax  =  <$  =  i  =  - v — --  meters 

x  B  4v  sin  0., 


If  we  now  consider  c  =  vX  the  above  becomes 


(A/2) 


6a  2  sin  0, 


meters 


and  about  what  we  would  expect  for  a  lens  operating  at.  a  wavelength  with  a  numerica1 
aperture  of  sin  0m  [31.  The  differences  between  (12)  and  the  expressions  of  optics 
are  the  factor  2  because  we  are  dealing  with  a  two-way  path,  and  a  multiplier  which 
reflects  half-power  points  of  the  correlation  function  and  not  the  ^auare  response  we 
have  assumed  here.  More  familiar  in  the  radar  arena  is  the  small-anqle  approximation 
of  (12) 


(d) 


meters 


Both  of  the  above  forms  should  be  quite  familiar,  but  here  they  have  been  developed 
assuming  nothing  more  than  a  monochromatic  interferometer,  a  particular  data 
collection  format,  and  matched  filter  processing.  The  interferometer  and  the  ma^rhpd 
filter  comprise  the  only  physics  of  Synthetic-Aperture  Radar.  The  major  efforts  ni 
SAR  deal  with  engineerino  aspects,  such  as  those  related  to  processino  which  wil1  ho 
outl ined  in  the  fol lowing  sections.  A1 1  we  have  to  do  j_s  engineer  these  s imp! e 
physical  principles . 

Just  as  we  may  process  the  data  in  the  azimuth  direction  with  a  matched  fffer  w n 
may  consider  our  interferometer  to  have  a  mul tif requency  bandwidth  RT  and  matched 
filler  the  optical- path  difference  returns  at  the  time  delav  or  range  of  interest. 
These  returns  may  be  compressed  to  a  pulse  of  length 


V  V  /  V  2 
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At 


(14) 


corresponding  to  a  range  resolution 


6 


r 


CAT 

2 


meters 


(15) 


In  practice,  we  will  have  compressed  the  transmitted  pulse  T  into  a  signal  of  extent. 
At,  or  achieved  a  compression  of 


C 


omp 


(16) 


equal  to  the  time-bandwidth  product  of  the  transmitted  signal.  The  returns  of  each 
PRF  line  must  be  compressed. 

The  item  of  note  in  this  section  is  that  the  processing  of  our  optical-  path 
difference  interferometer  data  can  be  accomplished  using  a  matched  filter  if  signal 
detection  and  relative  signal  strength  data  are  the  only  requirements.  The  radar's 
performance  is  determined  by  the  temporal  bandwidth  of  the  transmitted  waveform  and 
the  spatial  bandwidth  of  the  azimuth  or  aspect  data.  The  time/space  data  are 
compressed  from  their  measured  extents  by  their  time/space  bandwidth  products  during 
the  matched  filter  process.  If  the  frequencies  to  be  processed  are  centered  about 
zero,  as  in  (10),  we  are  in  a  minimal  sampling  requirement  format  and  find  the  need 
for  a  minimal  number  of  samples  equal  to  the  time/space  bandwidth  product  of  the 
signal  of  interest. 
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PHYSICAL  ARRAY  EMULATOR 

The  preceding  discussion  deals  with  the  after-the-fact  processing  of  data 
collected  by  a  coherent  measurement  system.  However,  the  processes  have  an  analog  in 
the  form  of  a  real  array,  as  a  real  array  can  be  built  which  will  yield  the  same 
results.  In  this  case  the  physical  location  of  each  element  would  be  known.  It 
need  not  be  on  a  plane,  but  we  will  assume  the  elements  don't  block  one  another. 

Considering  the  array  as  a  whole,  we  want  the  energy  from  all  elements  of  the 
array  to  arrive  at  a  resolution  element  in-phase  (constructive  interference  of  the 
multibeam  interferometer) .  Each  element  will  transmit  an  identical  waveform  which 
provides  the  desired  range  resolution.  The  desired  overlapping  is  accomplished  by 
providing  the  desired  waveform  at  each  element  and  then  adjusting  its  phase  and/or 
time  delay,  relative  to  a  local  reference,  to  offset  the  path  delay 
kr  (modulo  -  2n)  or  relative  time  delay  to  a  given  resolution  element.  This  energy 
is  scattered  from  the  target  and  returned  to  the  individual  elements  of  the  array. 
Upon  reception,  the  sional  at  each  element  must  again  be  adjusted  in  phase  or  time 
before  it  will  constructi vely  interfere  with  the  returns  collected  by  all  other 
elements  and  provide  us  with  the  full  array  performance.  To  probe  some  other 
resolution  element  we  appropriately  readjust  our  phasing  and/or  timing  at  each 
element  and  repeat  the  process  and  continue  in  this  manner  until  we  have  interrogated 
the  space  of  interest. 

In  the  SAR  we  have  removed  the  time  coincidence  restriction  and  have  combined  the 
two  phase/time  operations  into  one.  Other  than  that  it  should  be  apparent  that  the 
brute-force  matched-filter  SAR  is  an  analog  of  a  phased  array  with  its  elements 
positioned  at  the  measurement  points  of  the  SAR.  Both  are  interferometers.  The  real 
array  is  a  real-time  multibeam  interferometer  and  the  SAR,  indirectly,  is  a  multibeam 
interferometer  by  virtue  of  its  coherent  standard  for  all  measurements. 
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PROCESSING  CONSIDERATIONS 

The  preceding  discussion  deals  with  a  single  point  in  the  field  of  view  of  t h < 
radar  antenna.  In  practice,  we  are  interested  in  all  potential  target  positions  ir  - 
neighborhood  of  P  as  illustrated  in  Figure  2.  We  can  reasonably  assume  that  we  can 
adjust  our  antenna  pattern  so  as  to  view  all  points  in  the  area  X  (ranoo  and 
azimuth)  as  the  synthetic-aperture  length,  L  ,  is  traversed.  Thus,  our 

O  jr  I  I 

interferometer  will  collect  data  on  all  points  in  the  target  field  of  interest  durino 
its  travel  through  the  distance  L  p.  We  will  now  process  these  signals  to  detormine 
the  cross  sections  of  a  target  at  each  of  the  positions  in  the  neiahborhood  of 
interest.  To  avoid  considering  all  targets,  we  assume  a  linear  system  and  consider 
only  a  single  target  which  is  arbitrarily  positioned  in  the  field  of  interest. 

We  begin  by  assuming  a  brute-force  processing  of  the  data  with  a  range 
compression  factor  T  By.  Thus,  we  will  require  at  least  T  BT  complex  operations 
(multiply  plus  add  for  the  minimum  of  TB  samples)  to  compress  the  range  data  for  a 
given  resolution  element  for  each  PRF  line.  We  require  a  minimum  of  LB^  PRF  lines  or 
a  minimum  of 


(TB,)  (L  B.) 
T'  syn  A 


(17 


operations  just  to  process  the  range  component  of  a  given  resolution  element.  In 
addition,  we  require  a  minimum  of  ( L B ^ )  operations  on  the  data  after  it  has  beer 
compressed  in  range  to  accomplish  the  azimuth  compression  for  a  given  resolution 
element.  Considering  that  we  have 


resolution  elements  in  the  field  of  view  of  interest  we  arrive  at  the  minimal  number 
of  operations  per  channel  required  to  produce  the  desired  image, 

'LB,  •  'TB.)  (LB,) )  “r  ,-a 

r  a 

Tf  we  assume  TR^  1  and  use  Ill)  and  (13),  the  abovp  mav  bo  api  *  ex  imated 


11 
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(;> 


required  operations  per  channel  (range  processing  clearly  dominates).  The  numbers 
became  rather  large  for  practical  radar  values,  and  has  been  fought  for  years.  If, 
in  addition,  we  impose  the  continuous  strip  map  condition 


D 


a 


>  L 

—  syn 


we  arrive  at  an  approximate  requirement  for 

(%)*  <tv(£) 

operations  per  channel  per  image. 

Consider  the  following  radar  example: 

5 

R  =  0.5  x  10  meters 

A  =  3  cm 

6  =  1  m 

r 

5  =  1  . 5  m 

a  4 
D  =10  m 

IK  =  2  '?7 

i 

Ke  would  require  approximately  10^  operations  per  channel  to  produce  an  image  of  an 
area  of  (]0^)  x  (10^)  meters  squared.  The  minimal  synthetic-  aperture  length  would 
h e  1000  meters  or  about  5  seconds  for  a  platform  velocity  of  200  m/sec  (  ~  0.7M). 

10 

rhus,  we  see  a  basic  minimal  single-channel  processing  requirement  in  excess  of  10 
operations  per  second  for  a  real-time  strip  map. 

Processors  for  synfhet.ic  aperture  radar  data  have  evolved  alonn  fhe  line  nf 
■  qijf  f  i  a  1  ranrje  and  azimuth  processinn.  Part  of  this  may  be  optically  driven. 
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However,  we  can  see  a  rationale  for  this  by  considering  the  matched  filter  for  the 
point  P'  illustrated  in  Figure  2.  The  point  is  assumed  to  he  in  the  field  of 
interest  and  its  ideal  matched  filter  would  have  the  form 


-  i  '  ) 


i  2kr ' 


for  the  siqnal  form 


A  ( r  '  )  e 


■i  2kr  ‘ 


collected  along  the  data  collection  path.  Lookinq  at  the  filter  form,  we  note  that 
mathematically  we  may  represent  r'  as 

'r  3xt  |  t  Dy.  |yt  '  ’  '  " 

J  r  1 

0  c 

whore  r0  is  the  distance  to  a  reference  point  in  the  radar's  field  of  interest.  If 
the  higher  order  terms  can  be  ignored,  the  matched  filter  for  off-reference-pos i t.ion 
point'  ip  the  radars  field  car  be  viewed  as  linear,  target  position  dependent,  and 
perturbations  from  the  matched  filter  for  the  reference  point.  What  this  implies  i 
that  fhe  processing  can  proceed  by  first  filtering  all  da+a  along  L  usinq  the 
matched  filler  for  the  reference  pnin+.  The  residual  dafa  car  fhen  he  processed  ir 
some  maonnr  to  sort  out  the  desired  information  for  the  other  points  in  the 
np i qhhn>"hood  of  the  referencp  point.  We  will  look  at  some  historical  implemor4'- ‘ire 
hnluw,  but  le+,s  first  look  at  some  general  limitations  of  the  implied  ranqe-az imu? n 
~rocess inn . 


i  r.PSider-nq  range  and  cross-range  coordinates,  we  see  a  '  ’('nal  from  a 
M  aspert,  as  in  Figure  3,  would  have  a  matched  fi'tpr  form  app»e*  iina  *  ■  ■ 
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r'  (-  °H);  r  +  |^-  |xt'  f  r;r'  V  .  ...  ( 

3xfl  ' V  L  ’  'I 

From  +  ?i|  the  form  would  approximate 

r'  (0„)  :  r  +  |r'  |  y  .  ,jT  j  yt'  ♦  ... 

"xt" 1  r  '^t  r 

As  pointed  out  by  Brown  and  Fredricks  [”4]  the  residue  phase  can  be  processed  in  a 
range-azimuth  ( rectangul ar)  sequence  so  long  as  the  x^'  and  / 
representing  a  target  position  deviates  by  no  more  than  one-half  a  resolution 
element  as  viewed  from  various  aspects.  Likewise  foryt'  and  y^11  .  The  most  chanae 
occurs  at.  the  edges  of  the  field,  and  for  a  given  A6  we  can  impose  the 
f ield-of-view  limitations 


D 

~  A0  <  6  ( 

2  a 


(  t)A9  <  6r 

Using  (13),  the  above  become  the  expressions  of  Brown  and  Fredricks  f4], 


( 


1  • 


and 


5  6 


(32} 
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Returning  to  our  consideration  of  sequential  range-azimuth  processing,  i.e.  ,  we 
first  process  the  returns  in  range,  and  array  the  data  into  ranae  bins  for  subsequent 
azimuth  processing.  This  involves  compressing  the  returns  with  a  matched  filter  for 
each  range  bin  or 

D 


operations  per  PRF  for  a  range  swath  of  Dr  .  This  ranae  processing  may  also  be 
accomplished  in  an  analog  fashion  with  a  dispersive  delay  line  operating  on  chirp 
signal,  a  tapped  delay  line  for  other  signal  forms,  etc.  Under  this  condition  we  mav 
treat  Dr  as  cr  and  (31)  becomes 
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and  a  not  too  stringent  requirement  for  reasonable  resolutions.  However,  if  we  look 
at  (32)  and  impose  the  strip  nap  requirement  L  we  arrive  at 


nr 


6  5 

r  a 
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A  >  *2r 
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(35) 


(36) 


Th i s  latter  condition  can  easily  impose  resolution  values  in  the  meter  range  for  the 
assumed  sequential  ranqe-azinuth  scheme  of  processing.  We  can  contrast  this  to  an 
’ d°a 1  point-by-pnint  matched  filter  where  we  may  only  limit  the  half  angle  to  perhaps 
3b"  (because  range  and  azimuth  data  provide  the  same  resolution  characteristic  at 
90")  and  can  obtain  limiting  resolutions  on  the  order  of  *  .  We  need  not  concern 
ourselves  with  the  whys  of  the  above  limitations  as  we  will  later  look  at  procedures 
to  overcome  them.  They  can  be  explained  in  a  number  of  ways.  For  now  let  us  look  at 
now  processing  has  evolved. 


20 
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Radar  data  was  initially  processed  optically  by  recordinq  data  on  film  in  a 
sea -nr1  ranqe-a? imuth  format.  Optical  processlno  allows  compression  of  wideband  chir; 
squeals  with  good  range  resolution  and  with  azimuth  resolutions  on  the  order  of  those 
indicated  in  the  above  illustration.  The  most  successful  processings  have  been  those 
made  with  the  data  formatted  for  the  tilted  plane  processor  f5l.  While  optical 
pre^essing  has  many  nuances  let  us  leave  it  at  this  point  and  go  back  to  our 
description  of  the  data  form  and  the  digital  processing  impl i cat i ons .  We  will  return 
*  n  a  discussion  optical  processing  later  when  we  make  use  of  the  ?-D  Fr,uri«r- 
■  ransform  property  of  a  lens. 

I  e*  us  consider  herp  the  data  contained  in  a  range  bin  of  extent  •  as 
illustrated  in  Figure  4.  We  simultaneously  collect  the  data  for  all  target  pnsdmm 
in  this  raeae  bin  along  the  synthetic-aperture  length  l  . 

'he  matched  filter  for  an  arbitrary  point  in  this  limited  field  is 

ei2kr’ 
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Rather  than  imp! ementing  this  filter  for  each  Xf  let  us  look  at  the  phase  residue 
after  the  matched  filter  function  for  the  point  =  0  is  removed  from  all  data  in 
^he  range  bin.  This  operation  leaves  a  residue 


xx 


At  6 


12  k  (-/) 


(4 


tn  he  processed.  Note  that  samplina  considerations  have  not  been  explicitly 
addressed,  but  we  assume  sampling  is  adequate.  The  PRF  (azimuth  sampling)  of  the 
radar  may  be  averted  by  the  particular  hardware  implementation  selected,  as  will 
things  'l'ke  presuminq  during  data  processing.  Here  we  ore  addressing  the  minimal 
processing  requirements  of  some  optimal  data  collecting  implementation. 


If  we  1ock  at  (41)  we  see  the  residual  signal  retains  ?  bandwidth  of 


/  j_  | 

\  =  2/ls<l"\  =  v 

radians 

( 4? 

\  2"'  lx  1 

/  l  A  r  J  6 

m 

'  lx 

/  \  o  /  a 

f or  X  =1/2.  This  implies  (L  /<$  )  samples  and  processing  operations  for 

max  svn  F  v  svn  a  K  K  * 

each  of  the  (l  /5  )  azimuth  resolution  elements  in  the  range  bin  of  interest. 

These  operations  come  on  top  of  the  operations  required  to  remove  the  reference  point 
matched  filter  from  the  data.  The  reference  point  variation  can  be  removed  in  the 
a^aloq  portion  of  the  radar.  Digital  removal  of  the  central  matched  filter  imp1 ins 
■•continuous  strip  format) 


!  \ 


operations  to  process  the  azimuth  data  in  each  range  cell.  These  azimuth  operations 
are  in  addition  tn  the 
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operations  required  to  assign  the  data  to  appropriate  range  bin.  Tne  minimal  number 
of  operations  for  our  channel  strip  map  for  this  case  being 
(Operations/range  bin)  x  (#  of  range  bins) 


(1 


+•  L_  +  TBt) 
5a 


the  form  is  different  from  that  for  the  brute-force  case  as  are  the  number  of 
operations.  For  our  radar  example  given  above,  L?  /*  wi 1 1  dominate  TBt,  and  the 
range-processing  requirements  become  insignificant  relative  *o  azimuth  processinq, 
with  the  result  that  we  now  minimally  require  approximately 

(i 

operations.  Vie  have  actually  increased  the  required  number  of  azimuth-related 
operation,  but  have  reduced  the  required  number  of  range  operations  by  treating 
CL/ 5  )  azimuth  elements  as  though  they  have  the  same  range  filter,  i.e.,  always  lie 
in  a  single  range  bin.  The  original  brute-force  example,  a  direct  implementation  of 
the  mathematical  matched  filter,  assumed  each  resolution  element  required  its 
unique  range  matched  filter.  With  the  azimuth  processinq  a  dominate  factor,  we  can 
conceptually  explore  reducing  this  load  more. 

If  we  return  to  Mil,  we  see  that  for  small  9  (  r  a  weak  function  of  x)  the 
residue  function  for  a  signal  in  a  given  range  bin  takes  the  approximate  form 


where  :j  is  the  broadside  ranae  to  the  target  in  Figure  4  It  may  be  a  mean  target 
*'arge  for  other  squint  angles.  The  residual  siqna1  approx mid  +  es  a  constant  spatial 
t  require y  with  the  frequency  proportional  to  the  target  position  X  ^ .  The  desired 
processing  is  a  spectral  analysis  which  can  be  performed  batrhwise  using  an  FFT 
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itl  .j,.r  :  •  iia.  1  ho  (TT  w  i  1 1  bo  Addressed  in  wre  del.'!  1 1  be :  ■  w .  '  ;  >:•  pe.:  mm-.*  j 

trie  c  .  j i  »'o  !  number'  -if  : > I'i., i  b  operations  to  the  <>t  lc  •  •  f 


7  he  FFT  >'npl  '*">  ;L/-  M  sr  integer  power  of  2,  and  is  some  value  so  tha*  i 
*hc.  =  z  ■■  mu r h  swath  width  0„ .  The  number  of  azimuth  operations 


o  i  v  e  n  by  M3)  would  be  the  required  number  for  a  brute-form  matched  filter  approach 
-nurier-'t-ransform  operation  requires  a  like  numh°r  it  an  FFT  algorithm  is  not 
emp1oved.  Processed  w’th  a  FFT  algorithm  the  required  per  channel  operations  new 
■i  ppmx  ’  ma  he  s 


and  ranqe  compression  may  dominate  once  more. 

The  above  illustrated  an  efficiency  offered  the  azimuth  processing  when  a  common 
tag tor  is  removed  from  a  11  da*a  in  a  range  bin  leaving  an  azimuth  data  residue  which 
can  be  processed  with  an  efficient  algorithm.  A  similar  efficiency  tor  fhe  ranqe 
da +  a  may  he  postulated  it  the  ranqe  data  were  in  a  form  which  could  be  processed  by 
spectral  analysis.  This  is  the  data  form  produced  when  the  retch”  technique 
described  hv  ''a pat’  (6"!  is  employed. 

Thp  - h ’ r pod  return  f rnn  a  target  point  at  range  r* 

l?  v  -  2r'  )  ‘  .  ft  -  h'",2) 

"f  "  0  c  2  c 

■'  n i  iic  w’tn  a  re*r>rrre  '•iqndi  form 
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i  2::  [j  (t  -  2r)  t  ,  (t  -  2L)1) 
0  c  2  c 


rather  than  the  fixed  freouency  reference  implied  in  all  previous  illustrations.  rhe 
result  is  a  monotone  video  signal  whose  frequency  is  dependent  on  the  differential 
ranqe  (r'-r).  This  signal  has  the  parametric  frequency  form 


The  transmitted  BT  is  T  ard  a  resolution  c/2By  ^an  he  produced.  The  processing  of 
these  signals  is  a  spectral  analysis  and  an  FFT  may  he  employed,  i/e  arp  now 
considering  the  data  from  more  than  one  range  bin  and  the  range  swa4'h  imposed 
limitation  (31)  may  come  back  into  play.  Brown  and  Fredricks  were  workinq  with  such 
chirped  siqnals  and  the  rotary  table  SAR  emulator  at  the  Willow  Run  Laboratories  (now 
"RIM)  when  their  paper,  referenced  above,  was  published. 


While  we  may,  on  the  surface,  now  anticipate  more  reduction  in  our  processing 
!cad,  the  range  data  structured  as  above  and  implemented  with  an  FFT  requires  some 
tradeoffs.  We  recall  each  PRF  required 


brute  force  range  processing  operations  to  rompress  the  data  and  assign  it  to  one  of 
./■'  )  range  bins.  We  operated  under  a  swath-width  constraint  with  the  form 


2 

a 


and  were  quite  happy.  Processing  this  modified  range  data  wifh  a  :  ■  r . , 1  - 1 . .  r  . 
Fourier-transfcrm  will  require  the  same  number  of  operations.  However,  we  would  ’ike 
to  exploit  the  FFT  algorithm's  efficiency  in  our  spectral  analysis.  This  romes  at  a 
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T'v  TP  product  of  the  input  signal  dictates  the  minimal  size  7  our  f-R.  Ir  the 

example  above  we  assumed  30  for  a  TB  product.  He  could  implement  A  3p-p<  mi  (:(  r  ,t  r  r, 

cad  with  .-pros.  The  processing  would  separate  data  into  3?  one-meter  range  hi  ns. 

'Hi  3 -’-mpi'er  swath  applied  to  (311  does  not  impact  our  desired  l.su  'inuth 

resol  u* i  ur .  Th.e  prot^er1  is  tha*  we  have  been  forced  to  cons  trier  or  i  v  a  ne.--MW  nmjr 

swath  wher  low  tb  products  are  employed.  For  our  model  radar  wo  would  have  -n 
d 

It'  "TP  ;eparafp  L.O.'s,  with  each  representing  an  appropriate  range  del:,. 

Ho  wove-*,  let’s  proceed  and  note  that  we  now  require 

TBr  log,  (  7R  T  ; 


-ress  ’  ng  operat’  -Ts  for  each  of  subswaths  of  width  '■  TP.,1 


For  part  f-F!' 


■og. 


-  ':rr?,  +  'ors .  A  rather  curious  result,  says  that  the  processing  load  goes  down  as  we 
decrease  the  TR  product.  We  do  indeed  reduce  the  processing  burden  to  the  1 i m i t  of  a 
shor-  pulse,  IJ3~)  -  !.  However,  we  have  shifted  the  burden  to  an  extremely  complex 
server  *-haf  orobablv  could  not  be  implemented. 

!  *  we  take  *  he  *ack  r.f  minimal  L.O.'s  in  the  above,  we  find  (31  )  limits  pm  a 
•  n*  swa  * h  .noth  give  Pv  'example  numbers  assumed) 


u ■ 


w  id  f  h  : 


I- 


n 


•s 


/v 


5S 

•»  4 


(  .  arp  i  f.  va1  ues  ’his  would  be  about  1  #37  meters.  let  ns  select  TB  •-  )  ?H ,  arc; 
'<  ’  ■■auen* '  y ,  a  Hibswath  r  f  138  meters  for  the  assumed  range  resolution.  We  would 
’•p'.'jino  ,>i-ci,t  P,0  local  e « c i  1  latprs  and 


'V'Vv  v 
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log,  '12g)  -  7  "r 

r  Jr 

ranoe  processing  operations .  This  number  is  not  a  lor  dif+erent  'nr  +np  bnj'c-'i  rr 
'70  (r  /  }  indicated  fnr  the  brute-force  processing  of  an  arbitrary  ~PT  ■  TO 

si  goal,  assumed  above.  This  small  improvement  comes  at  the  expense  of  a  severe 

implementation  penalty. 

The  above  would  seem  to  indicate  that  for  moderate  reso1  ut  i  nr- , ,  the  sepa^'cr  r 
data  into  ranee  h’ns  for  subsequent  azimuth  processing  can  a**"ord  s ’cr i ' ’ cant 
process  -'00  efficiencies  over  the  brute-force  approach  to  the  ideal  matched  fi’ter 
operation.  The  azimuth  signals  in  a  given  range  bin  can,  in  general,  be  additional1 
formatted  and  batch-processed  using  an  efficient  algorithm  to  affect  PVPn  greater 
overall  processino  efficiency.  However,  a  formatting  of  the  range  data  into  a 
similar  form  for  batch  promssing  provided,  at  best,  mixed  results.  We  did  not 
nbsprve  these  mixed  results  in  the  azimuth  processino  since  they  were  not 
anticipated,  and  the  space  bandwidth  product  processed  was  only  that,  required  to 
prrduce  the  desired  map  rate,  i.e.,  a  real-time  strip  map.  While  we  may  consider 
breakina  the  azimuth  into  small  batches,  Appendix  A  shows  that  this  does  not  reduce 
the  digital  processing  load.  There  is  a  large  overhead  in  the  formatting  operation 
which  was  not  visible  in  the  range  example  given  above  as  the  complexity  was  forced 
bark  into  the  hardware. 


'•‘'•V  •  !  V,  I'f  ••■AHAR  :r? 


forma*  finq  of  »Mda»*  data  was  conceived  *• ,  H  •  '  ->  •  e*  p.-  ■  * 

'  i  *  ''-'Os  dps  by  R  rr  wn  and  Fredricks  f  4  1  dur :  "g  »hpir  »•*  *  •  w,*>'  •>»  ' 

*  -*■'  .  ;  1  h *-'*•' vr,  "AP  pmu’ato^.  Tr'  *  n  *  s  s.  r,  w.  -  •  •  •  •  •  •  ♦ 

••  •  >  i  i  T  is  for  the  concept ,  and  re*  a  [-irt.v.  e  :,<r  -i  \  nr<  ?.••*  e 

•  • ■.>  r.p  s  pr  va T  i  c  *  Wa  1  ker  arp  quite  p  :  egart  and  rr  o  phvs'r-*'  •  *  •  tun 

;•  -n.  rod  w'*f  *  !.<■•  ’f'  +  rpnuct  ion  nf  the  chirp  waveform  and  '■<  1  .•>  -  r  r<c::,i 

v  .  ,‘-a  ops/'  'jd’np  nf  the  concept.  The  ch’rp  w^vo^orr"  prnv: '*.(•'  a  - nr 
■  .<'.»>  ipp'v'ach .  but  is  not  basic  to  the  concept .  h  add  it :  •  w  ■ 

■  ■ r.o’n*  ,v'  >  •  So  :  -nv  :ded  later  in  ♦‘his  section. 

’ r  oc-v e’,:p  the  ba<  * '  c*  ‘•he  polar  format,  consider  the  target  rpfpron.  eg  dat 
*  ->r  cenr-ptrv  i '  ■  us  * '"afod  in  Figure  5.  We  assume  that  ♦argots  n*  ;n<fTs‘ 
*  e  rpH  m  *-np  neighborhood  of  the  origin  of  an  inertial  1  v  defined  f  / ,  tv4  i  p 1 
u- d  da* a  ,  r;l  1  oc*  i ,-,’i  occurs  from  a  remote  inertially  determined  poir.t  P,  w*‘ 
ally  np+  in  the  ■  x ,  v ;  plane.  The  observation  point  P  and  *he  z  axis  define 

f  h<’* ;  i-  a  i  ■  Ip,:  t  ion  piano  which  would  intersect  the  y  )  plane  as  shown. 

(,f*  makes  ,:r  angle  ’•  with  ♦he  y  axis.  The  1  ine-of -sight  distance  from  ih 

■  aatior  p.;  int  P  to  ♦he  reference  point  is  denoted  as  r  ^  ,  and  intersects  the 
./  plane  a ‘  an  i nstantaneous  observation  depression  angle:  .  Thp  1 ine-of- 
<■‘•<1  f  *  e.m  P  *n  ar-  arbitrary  target,  position  (xt,y,l  1  s  denoted  bv  r  . 


.<•  a  pe 1  oc :  ye  crr  woat  me-,  seem  a  weird  and  unconventional  selection  "f  >  n  <- 
7i *  •  - ■ s  ,  i . e . ,  a  Vf *-handed  basis  system.  However,  there  is  a  method  •  i 

i:r,r,,  no ro  cv •  do,it  as  wp  discuss  imaqe  forms  be  1  ow .  For  now  wr  w i 1  ]  , 

■  ■  -  * ;  r  ’  4 ;  *  c  in-,  w  •  ‘  h  nominal  movement  in  the  y  direction,  i .  e . ,  antenna 
>-v,,rri  s.  i  dn  rf  c  '•  1  lection  vehicle.  Using  the  coordinates  shown ,  the  imago 
,r ed  wit’  pr-  ' ■  n e  wh i ch  would  correspond  to  what  an  observer  would  see  1  r/  i 
>  : -  r  a  n,-;  c  '•  d»  nt  r.,',1 1  or*  ion  platform.  On  a  display,  range  may  bn  i  ‘‘v 

-  a ;  T : f  i  ,-,n  ,  with  far  ranges  at  the  top  of  the  d>spi  ov,  and  ♦ho  r-< .  s  ’  t  i  vr 

■  ,  r  on  *•.<  observer 1  s  If  ft..  In  a  strip  man  mndf1,  *he  imago  would  move 
fry.  *•>,.  d-'splav  vs  * i me  proceeds .  In  the  same  manner ,  the  ^ada r  so,.'  r 
v  ar.  o t' serf  >  on  f  he  re!  lection  pla  +  form. 


vV 


Data  Collection  Path 


lo 1 1 ection  Geometry  tjsei  U. r  Spot! i gh t 


30 


p7p  p  'ocjdnce  and  power  of  polar  f  f > r; '- .  f  •  o,  .  ;--p  ■’ 

1  if>r*  -  or  svstpr  os  i  nterfc»'onet>  r  w  ■ '  f  f  rr  , 

■~o  prop*,  ft  i  t'S  of  ,ln  ’  ♦e^f  prompt  0^  ’  S  'r,‘  .  ,i  K  i  . 

.  orient  i  l  c  perhaps  f  apar  •  <  *  ■■  i, 

.  -.H  ir-  ,i'4:c  •  i; r  1  •'  nea  r  s  V c  ’ :  oi,  •  •  ♦  ■  < 

■■  *•  7  .  s«*  i"‘ »  c  r.f  nonr*f  n*-s  ,  w :  »  '•  ■  v‘  .  ♦ 

•  '(•  i  '*•  •  •’(.•;<=  1 -t'l-noisr  .  rhr  »•.  • 

'  '  •  ,  ,1  .  •  •  I  ,  *  r>smi  t  *Ctj  r-r) rif .  I  .  ,  r  :  •  .  .-i 

■  *  1  *  •  •  •  ■  1  ’  ' : •  ,  f  h  p  v-p  t  u  r < .  1  ■  ■  •  ■  t~ 1  ■ 

'  .  ’  o'  '  ,  '■  ;  ’  (fd  u'-’W;  ‘  K  ’ 


w--  w; '  ’  receive-  a  '  f.t »k 


<_ 

p., os  an.  doing  hero.  This  ••  <:  -onpar, 


1  s  or  optical -path  o i *• r;-p 


1  'll  ‘  ifTOn.  r  hiifvcff.- 
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viewing  aspect,  i.e.,  different  9  ,  and  repeat  our  measurements.  We  have  not  required 
continuous  motion  nor  a  chirp.  The  only  requirements  are  a  knowledge  of  r  ,  ,  an d, 

as  we  will  illustrate  shortly,  the  angle  9  associated  with  eac'a  observation. 

The  rp^prrrrc  process  outlined  above  is  the  removal  of  the  optical-path  data 
representative  of  the  point  r  from  the  composite  of  data  collected.  We  have 
•~emoved  the  matched-  fi  1  ter  data  fo_r  point  (0,0)  from  al  1  data .  For  other  scatterers 
in  field  cf  interest,  we  are  left  with  a  phase  residue  which  is  dependent  on  , 

,  .  ,  :<r  and  y  Le*  us  attempt  to  characterize  this  residue. 

Expression  r  as  a  function  of  rn  and  the  various  parameters  we  find 

, ?  2  91/? 

r  =  .  r  sin  ,.iw  'r  cos  v  sin  0  -  ii<.  )  +(r  cos  c  cos  +  y.)1']  “  '53} 

0  O  t  0  /t 

Expanding  r  in  a  binomial  series  about  cq  we  find 
r  =  r  -  x  cos  c  sin  +  y  cos  cos  Q  +  ... 


and  in  the  paraxial  approximation  limit 


( r  -  r  ) 
0 


L-cos’l'S  i  n 


(55  ) 


As  once  stated  by  Walker  F8];  "The  form  begged  for  polar  representation . "  The  basic 
polar  parameters  being  and  6  with  the  -  coordinate  scaled  by  cos^  .  In 

this  illustration,  the  angle  coordinate  0  is  not  the  conventional  polar  coordinate 
a n r; i e  •  .  Wp  assume  that  the  data  space  cartesian  axes  correspond  to  those  of  the 
radar  space  as  illustrated  in  Figure  6.  The  data  from  any  measurement  made  at  angle 
"  is  perceptual  1 y  formatted  along  a  radial  line  in  this  data  space  at  a  radius  of 
■o'  .  .  The  radial  line  makes  the  same  angle  0  to  the  /  data  axis  as  the 
:'Glecfion  plane  intersection  made  to  the  y^  radar  space  axis.  The  data  space 
'artesian  npnrdi naf es  being  related  to  ,  v  and  by  the  t ransfermat irr 

f* v r  rPr,s  ions 


3? 


_  s 


( V  -  v  /  s .  *- ,  , 
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=  -  .  cos  p  sin  0  '6" 

and 

=  v  cos  p  cos  0 

We  assume  that  the  number  of  frequency  samples  along  a  data  ’ino,  such  as  that 
illustrated  in  Figure  6,  is  adequate  to  unambiguously  specify  the  phr.se  of  the 
various  range  cells  the  radar  space  as  viewed  from  that  perspective.  Moving  to  a 
new  perspective  will  provide  a  new  set  of  phases,  and  we  also  assume  ar  adequate 
sampling  in  the  0  direction  to  unambiguously  specify  the  change  betw^n  observation 
angles.  If  the  measurements  are  made  sequentially  from  a  moving  platform  it  should 

he  obvious  that  0  changes  and  the  data  line  will  not  be  a  radii  in  <.  ,  , 

space . 

Proceeding  with  our  frequency  and  angle  sampling  of  the  radar  space,  we  cmid 
eventually  (a  week,  a  second,  10  microseconds)  populate  the  data  space  in  the 

neighborhood  of  (  v  cos.,  8q  )  with  individual  target  data  samples  representati ve 

of  continuous  signal  with  an  approximate  form 

\  «'  T  [xd  xt  *  yd  V 

While  perhaps  rof  obvious,  this  formatting  to  first  order,  removes  the  data  wander  of 
figure  3  which  concerned  Brown  and  Fredricks.  However,  range  and  azimuth  have  los* 
their  distinct  relationship  to  data  axes.  To  extract  the  desired  cross-section 
information,  A  ,  we  may  want  to  perform  a  2-D  spectral  analysis  as  the  approximate 
signal  has  the  form  of  a  2-D  spatial  monotone.  Each  potential  scatter  position  is 
represented  by  a  unique  2-D  spatial  monotone.  An  efficient  spectral  analysis  which 
w’"M  simultaneously  extract  data  on  all  monotones  of  interest  is  a  2-D  inverse 
F>u>si<->r-trarsform  of  the  data  in  an  extended  aperture  of  the  da+a  space. 

Because  nf  the  form  of  a  2-D  inverse  Fourier  transform,  we  can  translate  our  data 
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a*is  refetence  to  the  position,  (  cos  (••  )  ),  and  incur  nr  mnro  ♦h,---  ■■ 

O  O  ’  0 

i material  constant  phase  offset.  The  processing  aperture  may  now  he  viewed  .>■,  *h: 
shown  in  Figure  7.  In  this  aperture  the  data  exists  as  samplings  of  the 
superposition  of  many  2-D  monotones.  The  samplings  having  nominally  occurred  a’r-'-- 
offset  center  radial  lines  such  as  those  indicated.  While  the  sampling  fnrv 
into  play  in  processing,  we  can  view  the  signals  as  continuous  and  use 
Fnurier-transf  orm  properties  to  provide  a  first-order  description  nf  +  hn  prone  - . 
image  characteri sti cs .  These  descriptions  will  comprise  the  remainder  of  this 
section.  Processing  details  will  be  considered  in  later  sections. 


We  have  assumed  a  linear  system  to  this  point  and  have  dealt  with  an  nrhi  +  mv  -  '  . 
positioned  target  in  some  neighborhood  of  radar  space.  We  have  defined  the  fire 
order  characteri sti cs  of  the  residue  signal  which  remains  after  its  formatting  in  fi 
particular  fashion  and  have  defined  a  processing  scheme.  Now  is  a  convenient  t-r, 
consider  we  do  not  have  one  signal  of  the  form  (69)  in  the  data  aperture  but  a 
composite  with  the  form 


Vi  e'  T  [xd  xtj  +  tj  J'tji 


We  car  process  this  signal  usina  the  2-D  inverse  Fourier  transform  with  the  resui 

S  'x.,  y.)  *  /°°  r  (2  A.  .  [x,  x  .  +  y  ,  y  .) 

o  i’  /i/  ■  tj  e  >  d  tj  J6  •/tj' 


i  4  TT 


-xd x,  *  yd  n1 


)  dxd  dyd 


In  the  above  x .  is  aligned  with  x^  and  y ^  with  y^  .  Interchanging  the  orders  of 
summation  and  intergration  we  auickly  arrive  at. 


d'J  1  x 


/ . 
i 


V 


1  ;J  t  >  < 


tj 


Vi’ytj 


The  signal  magnitude  at  an  image  position  (x^,y.)  being  proport ^ona 1  to  the  strength 
of  the  sinnal  Af  scattered  from  the  conjugate  geometric  position  in  radar  space, 
i . e , ,  we  have  an  "image"  or  map  of  the  scattered  field  amplitude. 
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In  practice  the  image  is  not  a  series  of  delta  functions,  but  has  a  resolution 
limited  by  a  finite  data  aperture.  This  finite  aperture  leads  to  an  image  positions 1 
unrerta inty 


Ax . 


=  Ax. 


ar„dT 


a  n  r| 


'Vs 


afydT 


where  Lfv4  1  and  1.1/  1  are  the  data  extents  in  the  orthogonal  xH  and  yd 

data  directions.  If  we  consider  the  broadside  case  we  have 

C 

"*yt  ’  2Bt  cos  ? 

I  o 


and 


‘"xt  2(v  cos  ip  A9) 
o  o' 

where  is  the  spread  of  frequencies  employed,  .  q  is  the  broadside  depression, 
is  the  center  frequency  and  A9  is  the  angular  extent  of  the  data  considered.  The  y+ 
resolution  is  (1/cos  ^  worse  than  the  C/2Bj  indicated  previously  for  an  ideal 
matched  filter  and  is  merely  a  reflection  of  the  slant-ranqe  to  ground-ranqe 
relationship  illustrated  in  Figure  8.  The  cross-range  resolution  is  also  scaled  t>or 
the  indicated  matched  filter  resolution  (>  /2 AO  )  by  the  same 

(1/cos.^  factor.  This  may  be  somewhat  confusing  until  we  consider  tha*  the  previous 
resolution,  (  •  /2Au  ),  involved  the  change  in  aspect  in  the  c  ol  1  e_c  t  i_on  snam.  Thr 
'.e  in  ( 7h  1  is  a  change  in  projected  space.  For  this  example,  the  collection  -lot  form 


P U •  II ' f 1  V«  I  «|  >.D >  I  i.i  >11  ■  » ■  i ■  •  m "j v  11 


W '  ',r  v  it  y.1  v  *  w  -  i  -  »  ' 
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would  have  only  viewed  the  target  through  an  aspect  of  cos  and  our  original 

matched  filter  guidelines  remain  inviolate.  The  data  collection  anqle  change  is  less 
than  the  formatted  for  any  .  unequal  to  zero. 

The  above  presented  only  a  digital  perspective  for  the  processing  of  the  data. 

We  have  conceptually  formatted  or  arrayed  the  data  for  a  large  0  chanoe  in  an 
annular-like  strip  such  as  that  illustrated  in  Figure  9a.  The  discussion  above 
centered  upon  the  processing  of  a  subaperture  of  this  data  with  the  processina  axes 
aligned  with  the  format  axes  as  shown  in  Figure  9b.  The  result  is  an  image  with  an 
(x,  y)  orientation  corresponding  to  the  broadside  (0  =  0°)  view  of  the  target.  a 
first  order  the  image  has  a  form  which  is  insensitive  to  the  actual  portion  of  the 
total  data  record  being  processed.  The  data  may  also  be  processed  against  any  other 
set  of  reference  axis,  for  example  as  those  of  Figure  9c  where  the  y^  axis  is  aliened 
along  the  0  direction.  In  this  case  the  image  will  have  a  perspective  of  the 
target  area  corresponding  to  the  view  from  the  nominal  angle  of  the  data  collection. 
The  choice  of  processing  axis  is  basically  not  process  dependent,  and  the  choice  may 
well  be  made  for  other  considerations.  Noncoherent  addition  of  images  to  reduce 
speckle  may  indicate  a  preference  for  the  nonrotating  image  produced  by  the 
orientation  of  Figure  9b.  An  image  providing  the  observer  a  target  perspective  would 
be  produced  by  the  orientation  of  Figure  9c. 

The  above  would  seem  to  indicate  that  we  have  eliminated  all  our  problems  and  can 
proceed  to  consider  processing  details.  However,  we  should  keep  in  mind  a  couple  of 
caveats.  The  foremost  is  the  actual  structuring  of  the  data  into  forms  that  are 
useful  in  our  conventional  transform  implementations.  This  aspect  will  be  addressed 
in  the  following  sections.  The  second  is  that  we  have  only  dealt  with  the  paraxial 
limit  of  the  binomial  expansion  of  r  about  r^  .  In  practical  terms  we  have 
considered  only  the  phase  residue  for  points  in  the  "near"  neighborhood  of  the 
reference  point  after  removal  of  the  matched  filter,  r  ,  from  all  data.  Factors 
which  can  perturb  the  neatness  of  the  above  are  contained  in  the  higher  order  phase 
terms  which  were  ignored.  These  factors  contribute  effects  analogous  to,  but  not 
one-to-one,  the  classical  aberrations  of  optics.  Distortion  is  analogous  and  is 
inherent  in  the  terms  of  (65)  which  were  ignored.  The  usual  chirp  spotlight 
implementation  also  introduces  a  distortion.  While  the  chirp  provides  a  convenient 
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means  of  scanning  a  bandwidth  in  a  known  and  regular  manner,  it  does  disrupt  nnr 
interferometry.  The  return  from  a  point  at  range  different  from  the  reference  poin* 
is  offset  from  the  reference  frequency. 

!  have,  in  fits  and  starts,  been  attempting  to  sort  out  these  effecfs  for  a 
number  of  years.  The  details  are  tedious  and  cumbersome,  and  as  noted,  the  forms  dn 
not  lend  themselves  to  classical  optical  interpretations.  However,  I  feel  I  am 
somewhat  near  and  this  report  may  serve  as  a  prod  to  go  back  and  attempt  to  generate 
some  sort  of  report  which  at  least  summarizes  activities  to  date,  and  af  best  may 
provide  definitive  insight  into  the  influence  of  the  ignored  terms.  Thp  forms 
considered  in  this  associated  analysis  will  also  allow  for  rapid  evaluation  of  fhe 
influence  of  deviations  from  the  ideal,  i.e.,  velocity  measurement  errors,  influence 
of  formatting  to  aircraft  angle  rather  than  ground  angle,  etc. 
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OPTICAL  PROCESSING 

To  this  point  we  have  dealt  with  the  general  collection  and  formatting  of  rada>- 
data  into  a  form  where  it  may  be  processed  to  produce  an  image  of  the  scatter  points 
of  interest.  The  primary  descriptions  of  processing  have  been  digital  based.  This 
section  is  included  to  provide  an  introduction  to  optical  processing.  The  polar 
formatted  signal  has  taken  a  form  which  may  be  processed  using  a  ?-D  inverse 
Fourier  transform.  If  the  data  are  physical ly  mapped  onto  a  coherent  ?-D  spatia' 
light  modulator  the  signal  phases  will  appear  as  ?-D  phase  variations  nr  a  nnmina’ly 
planar  coherent  optical  wavefront.  A  simple  lens  will  accomplish  a  2-D 
Fourier  transform  and  can  be  used  to  process  this  data  from  [9]. 

The  traditionally  optical  processing  of  radar  data  has  not  been  of  a 
Fourier-tansform  nature  and  a  much  more  complex  processor  is  involved.  The  spatial 
light  modulator  has  traditionally  been  photoqraphic  film  with  the  data  arrayed  in  a 
range-azimuth  format.  A  scaling  in  both  ranae  and  azimuth  is  required  to  record  the 
data  on  reasonable-sized  film  formats.  For  a  typical  azimuth  signal,  (38),  we  recall 
the  phasor  form  (37)  represents  a  real  video  signal 

At  cos  2  k  [R2  +  (v  -  xt)2]1/2 


Th i s  signal  is  amplified  and  used  to  expose  a  photographic  film  in  a  range  bin 
somewhere  across  the  films  width.  The  exposure  is  adjusted  so  that  the  processed 
film  will  have  a  transmission  coefficient  proportional  to  the  strength  of  the 
recording  signal,  (77).  If  this  processed  film  is  illuminated  with  a  coherent 
optical  phasefront,  of  sufficient  extent  to  encompass  the  data  collection  length,  the 
film  will  behave  as  an  optical  zone  plate.  The  focal  length  of  the  zone  plate  will 
be  proportional  to  R.  Distance  and  wavelenqth  scalinqs  both  come  into  play  to  scale 
the  zone  plate  focal  length  from  the  radar  *-ange  represented.  '  nr  r . , .  ,  ‘  : 1  r-  .• 
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plate  occurs  at  a  position  along  the  data  track  scaled  from  X^.  The  zone  plate  has 
one  dastardly  fault  in  that  both  a  real  and  virtual  image  are  produced.  A  virtual 
image  reconstructed  wavefront  appears  to  be  diverging  from  a  focal  position  behind 
the  data  film.  A  real-image  reconstructed  wavefront  appears  to  be  converging  to  a 
focus  in  the  real  space  in  front  of  the  film.  The  two  images  are  a  mani festat ion  of 
an  ambiguity  of  the  data  (77),  i.e.,  it  could  have  been  generated  as  we  postulated, 
or  in  some  other  experiment  where  the  sampled  energy  is  converging  to  the  target 
position.  In  the  optical  processor  this  ambiguity  is  handled  by  introducing  a 
spatial  carrier  which  forces  the  ambiguous  images  to  propagate  through  the  optical 
train  along  different  paths  with  respect  to  the  optical  axis  of  the  processor  and 
allows  the  use  of  stops  to  block  one  of  the  images.  The  video  of  a  chirp 
transmission  will  have  the  same  form  as  (77)  and  indeed  a  dispersed  chirped-  range 
pulse  can  be  recorded  on  film  and  processed  optically  because  of  the  same  zone  pla+e 
properties.  The  zone-plate  property  of  a  recorded  chirp  signal  explains  the 
traditional  use  of  chirp  waveforms  in  SAR  radars  based  upon  optical  processing. 

Other  waveforms  do  not  produce  smooth  optical  wavefronts  which  can  be  compressed 
optically  and  the  ranoe  data  must  be  compressed  before  it  is  recorded  on  film.  The 
aforementioned  tilted  plane  [5]  processor  uses  a  combination  of  a  one-to-one 
spherical  telescope  and  a  demagnifying  cylindrical  telescope  to  process  appropriate! y 
scaled  radar  space  data.  The  processing  is  the  2-D  focusing  of  severely  astigmatic 
s ignal s . 

The  Fouri er-t.ransform  (or  inverse-transform)  property  of  a  lens  can  be  explained 
in  terms  of  its  focal  properties  [10].  A  lens  will  focus  the  energy  in  a  planar 
wavefront,  incident  on  its  entrance  pupil,  to  a  diffraction-limited  spot  in  the  focal 
plane  of  the  lens.  The  focus  position  is  determined  bv  the  direction  cosines  of  the 
incident  wavefront.  For  a  monotone  such  as  (69)  a  linear  film  exposure  produces  a 
diffraction  grating.  An  incident  planar  wavefront  will  be  scattered  into  the  ?  1 
orders  of  this  grating.  Each  order  will  be  a  planar  wavefront  which  is  focused  to  a 
spot  by  the  lens.  Again  an  offset  is  required  to  separate  the  ambiguous  orders 
(imaqes).  In  digital  processing,  the  properties  of  quadrature  signals  are  used  to 
rpmove  thp  ambiguous  signal. 
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The  nuances  and  limitations  of  optical  procession  are  detailed  subierf  r  th 
own  right.  This  section  has  been  included  to  provide  an  1 ntrodur t ^ cm  to  The  run, 
of  optical  processing  and  not  to  debate  its  merit  for  radar  processing.  To  date, 
fl'lm  has  been  the  sole  media  for  recording  the  data  from  a  widr-ranqe  swath.  As 
spatial  light  modulators  improve,  optical  processing  of  radar  da+a  in  real  tniit  - 
become  practical.  However,  even  as  our  capability  improves  we  will  have  to  addrr 
topics  such  as:  the  practical  dynamic  range,  the  accuracy  of  this  analog  process 
and  the  realization  of  practical  apodizations  to  control  sidelobe  levels.  We  lea 
our  discussions  of  the  general  properties  of  optical  processing  aT  this  poind.  \! 
will  return  to  it  as  an  element  in  a  hybrid  back-pro iect ion  process  which  will  b>- 
described  later. 
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H I SCRETE  FOURIER  TRANSFORMS 


The  ever  increasing  availability  of  digital  processing  power  has  lead  to  its  nea 
exclusive  use  in  modern  SAR  implementations.  The  conceptual  polar  data  space  of  'he 
spotlight  radar  is  no  exception.  In  this  section  we  will  look  at  the  mathematic,)1 
nuances  of  processing  the  pol  ar-formatted  data  form  and  the  motivation  for  var-'ous 
implementation  forms  such  as  the  discrete  Fourier  transform,  and  its  efficient 
implementation  in  the  FFT  algorithm.  The  designation  "Fourier  transform"  is  used  i ” 
a  1posp  sense  with  preservation  of  the  observable  characteristics  ef  a  Fourier 
process  and  not  necessarily  the  precise  mathematics. 


Ue  recall  that  our  interest  is  with  the  extraction  of  information  from  da*,;  w ;  *  t 
the  <-psidue  form 


Mx,y, 


i  [w  X  +  w  y] 

k  x  y  J 


where  x  K  v  are  cartesian  coordinates  mappable  via  a  polar  transformat ion  to  the 
radar  frequency  and  a  position  in  radar  spacp.  The  Wx 1 s  and  Wy’s  in  a  compos -*0 
signal  are  to  first  order  proportional  to  distances  in  radar  space.  To  determine  fh 
strength  of  a  return  from  a  position  represented  by  (Wx,  Wv )  we  ;ie  r  fern  ,j  -  ■> »  .  •. 

:,er.;*’)n  on  'die  >-esj  j,,e  pvil  ,  i.e.. 


-i  [w. 


3  . 
1 


•V  V 


S  ( x ,y  )e 


+  w  y ) 


l.'e  recall  that  this  matched  filter  form  is  ideally  limited  to  a  small  neighborhood 
about  some  reference  position.  One  can  look  at  (79)  and  say  that's  a 
"Fourier  transform.”  Mathematical ly  it  is  not  quite,  as  we  see  a  minus  siqn  in  thr. 
k.e>ral  which  is  formally  associated  with  the  inversion  or  inverse  Fourier  transforms 
Also,  we  ideally  have  infinite  limits  of  integration  and  we  are  miss'nq  a  scale 
factor  of  1  /?  '  [11].  However,  all  this  is  a  "so  what'5"  It  is  the-  form  of  (:)'  fl.i 

:s  nf  interest.  The  result  will  have  all  the  properties  of  an  apodwed 
Fourier  transform,  savo  fnr  pprhaps  a  scale  factor  and/or  an  axis  invetxirm.  Ho- 
factor  is  not  important  as  all  tarqptc  are  treated  equally  and  rr  1  .*  *  i  ve  v  -  '  v 
w f 1 1  rema i r .  If  the  inversion  of  axis  is  a  problem,  wo  ran  wo>i-  with  the  n'h,  i 
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sideband  of  our  analoq  IF  chain  and  the  phase  2k  (r  -  r'L  It  is  not  important  ;t 

o 

(801  is  a  Fourier  Transform  nr  not!  We  will  illustrate  below  that  f!ii  ^hr> 
techniques  developed  for  Fourier  processing  are  applicable  to  179!  and  its  bas,r 
properties  are  the  same  as  those  of  a  Fourier  transform.  Let's  call  if  a 


Goinq  back  to  179),  and  for  purposes  of  illustration  considering  only  one 
dimension,  we  have  a  signal  form. 

SQ(w)  =  s .  (x)e  "iwx  d.x 

We  need  only  recall  e  "lwx  serves  as  a  good  approximation  to  the  ideally  mashed 
filter  for  one  of  the  many  possible  S.(x)  making  up  our  composite  signal.  The 
mathematical  descriptions  (79)  and  (80)  indicate  cumbersome  processes  in  which  fh‘- 
presence  of  each  possible  signal  has  to  be  probed  with  a  new  "kernal"  0  ‘,w*  .  !  >- 
us  seek  an  efficiency  in  a  digital  approximation  to  (80)  by  first  noting  tha*  f Sr'' 
is,  in  fact,  a  limiting  definition  of  a  summation  process  (Archimedes  is  credited 
with  the  conception) 

S  (w)  =  1  imit  ::  5.  (x)e  1  wx  I  Ax  . 

u  i  1  I  J 


In  this  fnrm  we  have  put  no  restriction  on  the  \x  .  1 S  other  than  they  all  tend  tn 

J 

zero  before  the  summation,  r  ,  is  equivalent,  to  our  integral,  /  .  Also,  there  is 
also  no  restriction  on  the  extent  of  signal.  Now  let's  sav  we  don't  ao  to  the  li”i 
of  (81 )  and  are  happy  with 

(Si  (x  )e  ’’;x '  1  ;.x  .  1 

x.  ] 

J 

where  N  reflects  some  finite  number  of  samples  over  the  field  of  the  signal.  A 
*irst  implementation  imposition  general1'/  comes  in  the  form  of  a  regular  sampling, 
wi^h  thp  sampling  rate  tied  to  the  frequency  rnntent  of  the  signal.  In  the  idml 
'  1 r”  i  *"  of  minimal  sampling,  we  assume  v.r  have  filtered  and  formatted  nun  si  grain  ,f 
absolute  bandwidth  B  in  frequency  space  into  the  zero  centered  hand 


Yu 


B 

2 


The  maximum  frequenrv  content  of  any  component  of  the  signal  is  B/2.  fThe 
sorting  of  negative  freguencies  requires  the  quadrative  channel . 1 
With  this  formattinq  we  can  sample  at  a  rate 

_  ] 

'sample  B 


and  satisfy  the  minimal  sampling  requirement  of  at  least  two  samples  per  wavelenr,h 
of  the  highest  frequency  component  in  the  signal.  At  this  sampling  rate  (82'  become 
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The  form  of  (851  has  become  a  little  more  structured  and  we  may  also  drop  the  scale 
factor.  The  minimal  number  of  samples  required  is  the  space-bandwidth  produc*  c(  the 
signal,  but  this  minimal  number  allowed  only  with  the  signals  in  the  band  of  (83'. 

The  above  is  beginning  to  introduce  some  structure  into  our  estimation  of  iw’1  , 
but  it  remains  a  tedious  process.  Let  us  now  see  whaf  we  may  accomplish  by 
introducing  a  structure  into  the  w  space.  For  one  thing  wo  know  that  i;'s  of  ir'i  ><-  * 
are  limited  by  input  bandwidth.  We  also  know  f  r>  <m  nir  matcher  -  f  i  1 1  or 
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e  in  -jncerta i n ty ,  and  any  input  frequency  may  only  be  determine  :  '  <■>  an 

e  ■-  •  lint. 


because  of  our  limited  data  aperture.  Stated  another  way,  for  the  form  of  (8^'  we 
cannot  expect  significant  changes  in  the  value  of  s  (w)  on  a  scale  less  than  the 
given  above.  This  being  the  case,  let  us  restrict  our  frequency  space,  w, 
interrogations  to  regular  intervals,  with  (88)  the  upper  bound  on  the  lenqth  of  these 
intervals.  For  our  signal  we  can  assume  w  =  0  is  one  position  we  wish  to  probe  arg 
may  choose 

“k  =  k;”  =  2'  U  :'9' 

with  the  result  (85)  may  be  manipulated  into  the  form 

N-l  -i  ki  dr/' 

S0(k/,w)  -  Z  S  (j.'.x)e  N 
j=0 

usina  various  of  the  relationships  expressed  above.  Inspection  of  )  shows  it  has 
a  periodicity  of  N  with  both  k  and  j.  These  are  manifestations  of  the  process  we  are 
performing. 

Let.  us  consider  the  k  periodicity.  The  signal  c  (k  •  w)  is  an  estimate  of  the 
composite  strength  of  frequency  components  in  a  subband  in  the  original  data  set. 

T h i c  frequency  component,  which  is  real  in  our  case,  could  also  have  been  the  msub 
o*  an  inadequate  sampling  of  a  higher  siona.l  frequency.  The  process  does  not  know  if 
the  sampling  is  adequate  or  not,  and  for  that  matter,  does  not  care,  and  will  provide 
an  output  for  all  ambiguous  signals  which  could  produce  the  sampled  data  set.  These 
ambiguous  possibilities  are  the  periodic  repeats  in  w  space.  We  know  we  have  a 
h.jr.rtl  ini  t-p/i  signal  and  as  such  only  nepg  consider  fj  intervals  r.f  extent 
i/[  in  a  rr(.quencv  band  of  oxtont  R.  *  he  process  sort  s  o;i*'  la  *  a  set  •  !  bin  iw  1  !  ‘  1 
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•i  ''ito  'i  bins  of  bandwidth  .tv.  With  this  sampling  of  data  and  frequency  space, 
we  see  a  need  for  a  minimum  of  N  samples  in  each  space. 

The  periodicity  in  j  can  be  explained  by  removing  the  summation  ir  y,  while 
retaining  that  in  W.  We  now  have  a  form 

.  2mxk 

-l  - - 

S0(k;.w)  =  /S.(x)e  L  d^  (91' 

Ir  this  form  S  (kv.w)  is  the  k^  coefficient  of  a  complex  Fourie--  spries  which  we 
can  use  to  represent  the  function  S.(x)  in  the  interval  L  [12],  i.e., 

.  2mkx 

S i ( x )  =  I  S0(kAw)e  L  >92) 

-oo 

Outside  the  interval  L  this  S^fx)  representation  has  a  period  of  L  or  N  X.  Thus, 
each  interval  L  contains  an  identical  signal,  but  again  we  don't  care  what  the  signal 
outside  the  interval  may  be  as  we  will  treat  that  in  some  other  processing,  if  at 
an  . 

With  this  insight  we  see  (93)  provides  us  with  estimates  of  the  spectral  consent 
of  a  signal,  S^fx),  at  discrete  points  within  the  frequency  band  of  +-he  signal. 

These  estimates  are  in  fact  the  coefficients  for  a  finite  Fourier  series  which  may  be 
used  to  represent  the  signal  in  this  interval.  We  have  a  limited  number  of  terms  in 
this  series,  hut  have  seen  mathematical  proof  that,  for  a  aiven  number  of  terns,  a 
finite  Fourier  series  provides  a  minimum  mean  square  estimate  to  S  .  (x ) .  However, 
we  can't  locate  a  reference  at  this  time. 

Cleaning  up  90)  a  little,  we  arrive  at  a  discrete  form 

N-l  2 n 

S  (k/.w)  =  T.  S .  ( jAx  )e  _1  ~n  jk  (93) 

j  =  0 
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which  provides  an  estimate  of  the  composite  signal  strength  in  some  suhhanri 
■  w  of  the  total  signal  bandwidth.  Expression  (93)  is  denoted  as  a  discrete 
Fourier-transforn.  The  form  d'scussed  above  centered  around  N  samples  in  frequ£in<", 
space.  Morn  samples  may  be  considered  with  more  processing.  The  minimal  sample  ’ 
N  reflects  the  restraint  physics  (bandwidth  and  sample  length)  places  on  +ho  ra+‘> 
change  of  S  (w) . 

The  complexity  of  (93)  can  be  illustrated  by  writing  the  N  possible  forr'  rf 
1 V  w  1  in  a  matrix  form 


s 

(i7 

s.  (1 ) 

0 

1 

S 

(2) 

c  /  *3  3 

0 

= 

3  ,  . 

kj 

= 

i 

• 

Jo 

(N) 

i 

S.  (N) 

l 

mm 

_ 

.•d'ere  a,.,;  is  an  M  x  N  matrix  with  elements  of  the  form 

kj 

kj  N 

The  array  indicates  the  need  for  N  complex  operations  for  each  of  the  N  possible 

p 

samples  in  frequency  or  on  the  order  of  N'’"  complex  computational  operations  to  rrriu 
*  he  data  set.  The  mathematics  of  Cooiey  and  Tukov  F 1 3  "1  and  the  efficiency  of  the  F 
algorithm  can  be  illustrated  by  considering  the  form  of  (94)  when  N  is  a  power  or  2 
In  this  case,  the  regularity  of  fa,  .]  and  its  modulo  ?  n  insensitivity  permit's  the 

kj  J 

di aaona 1 i zati on  of  the  matrix  in  1 ogpi  stages.  The  first  stage  of  this  process  i 
the  addition  of  row  k  to  row  k  +  N/?  and  the  subtraction  of  row  k  +  N /?  from  row  k. 
The  resulting  form  becomes 
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where  S  M )  s  are  permutations  of  S  s  and  S .  '  (  i )  s  are  now  sums  and 
o  o  i  ' 

differences  of  the  original  data  signals.  There  are  still  N  rows,  hut  we  now  have 

l 

two  arrays  ef  order  (N/2)  x  ( N / 2 1 .  Another  stage  of  operations  on  the  fa  kjl  and 

M 

ra  kj"1  residue  matrices  will  result  in  4  arrays  of  order  fN/4)  x  (N/4).  The  processes 
can  continue  to  iterate  through  log^N  such  stages,  until  only  diaqonal  elements 
remain  in  the  N  x  N  matrix  relating  N  outputs  to  their  unique  N  term  weiqhtinqs  of 

O 

the  input  data.  The  process  exploits  the  fact  that  the  N'  processes  indicated  in 
(94),  many  times  involve  cases  where  a  aiven  data  is  operated  on  by  the  same  weight- 
inas.  The  FFT  tends  to  collect  terms  involving  these  common  weightings  early  ip  the 
process  and  then  fans  out  the  result  into  latter  staqes  of  the  process.  An  N  point 
FFT  can  be  implemented  with  an  array  of  (N/2  log^  N)  Butterflies  (Figure  10).  A 
Butterfly  accepts  two  complex  input  signals  and  uses  one  complex  multiply  and  two 
complex  adds  to  deliver  two  complex  outputs  to  the  next  stage  of  the  matrix 
di agoni zat.ion  process.  The  total  process  being  completed  using  the  order  of 


N_  log-  N 
2 


operations  vs  the  tr  indicated  by  (94). 

The  efficiencies  of  the  FFT  algorithm  become  even  more  pronounced  when  a 
two-dimensional  filter  is  desired  and  we  have  to  consider  a  2-D  discrete  Fourier 


transform  form 


S  (kAo:,  mAm)  =  I  I  S.(jAx,  kAy)e 
o  j  1  1 


1  (kj  +  1m) 


2  2 

Here  a  brute-force  implementation  indicates  N  operations  for  each  of  the  N  possible 

permutations  of  (k,m)  where  each  has  the  range  N.  An  FFT  implementation  is  usually 

in  the  form  of  a  series  of  operations.  The  first  is  an  FFT  of  one  data  variable 

(,i  or  l!  for  each  of  the  other  variables.  This  is  the  analog  of  sorting  the  data 

into  range  bins,  and  may  indeed  be  a  range-binning  operation  for  a  given  radar 

2 

processing  operation.  This  still  leaves  N  input  data  points  which  are  then 
processed  with  a  second  series  nf  FFT's  with  the  roles  of  i  and  1  reversed  the 

azimuth  processing  of  range-binned  data).  Fach  sequence  requ’nrs 


✓  ,  S,  /  .  •’  . 

V.V  v. 
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2 


1og2  N 


Fi^terf  ly  operations,  or  a  total  of 

N2  log2  N 


for  the  complete  filter.  This  compares  to  the  operations  indicated  for  a  brute- 
force  implementation.  The  motivation  to  use  the  FFT  algorithm  is  obvious. 

Tn  this  section  we  attempted  to  provide  an  introduction  to  digital  SAR- 
matched-f i 1  ter  operations  which  take  the  form  of  a  Fourier  kernal,  and  aive  rise  to 
an  indicated  need  for  "Fourier-transf orm"  processing.  In  addition,  we  attempted  tn 
illustrate  the  order  of  complexity  of  a  digital  implementation  of  such  a  filter.  We 
illustrated  the  mathematical  basis  of  such  filters  and  efficiencies  which  can  accrue 
when  regularities  occur. 
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T NTFRPOLAT ION 

In  thp  preceding  section  we  observed  that  the  pol ar-formatted-SAR  siqnals  r.n'  be 
processed  digitally,  employing  a  discrete  Fourier-transform  algorithm.  Relative'1'/ 
e'ficient  FFT  concepts  may  also  be  employed.  However,  the  FFT  implementations  rcre 
with  the  burden  of  sample  lengths  and  more  importantly  the  need  fnr  regular 
samplings.  For  equal  resolutions  in  range  and  azimuth,  the  data  input  to  a  P-H 
FFT-based  process  has  to  be  in  the  form  of  data  values  at  points  on  a  regular  ?-0 
rartesian  grid.  The  data  collection  produced  values  which  in  genera1  ray  be  ranger1 
spaced,  but  are  nominally  along  obligue  straight  and  nonparallel  lines,  fpnterirg 
the  processing  aperture  on  some  segment  of  the  data  produces  the  reouirement.  <or 
estimates  of  the  data  values  at  points  on  a  cartesian  arid,  given  the  data  values 
along  lines  which  appear  to  diverge  from  a  common  remote  point.  We  assume  tha*  an 
oversample  solution  is  not  valid,  i.e.,  we  collect  so  much  data  that  each  cartes'/ r 
point  is  assured  of  a  data  point  in  an  acceptably  near  neighborhood.  We  could  thpr 
select  the  data  we  desire  and  throw  the  remainder  away.  Rather  than  bru+e-tnrcn, 
let's  just  assume  we  have  only  an  adequate  sampling  in  a  radial  line  format. 

One  method  for  the  estimation  of  grid-point  data  is  an  interpolation  based  upon 
the  data  form.  To  good  approximation  the  data  form  indicates  an  inverse 
Fourier-transform  processing  procedure.  Thus,  without  stretching  the  imagination  *o 
far,  we  car  say  the  data  itself  has  the  form  of  a  Fourier  transform  of  something 
which  closely  approximates  the  image  we  will  produce.  By  placing  a  few  restrictions 
on  the  data  sets  we  can  illustrate  an  interpolation  process. 

Consider  the  usual  frequency-time  Fourier-transform  pair  with  a  limited  time 
sample  transform 


F(.  ) 


1/2  i  t 

/  e  '  f(x)dx 
-T/2 


and  i*s  inverse 
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f(t) 


1 

/2tt 


-  lot 
e 

-  30 


F(^)doj 


We  have  produced  estimates  of  F(w)  at  a  series  of  points  we  may  designate  as  w., 
we  really  want  estimates  of  points  we  may  designate  as  ;•  ^ ,  i.e.,we  simply  go  bach 
and  find 


F 


/  T/2 
-T/2 


f(t)e 


V 


dt 


We  may  not  he  able  to  do  this  physically,  but  let's  do  it  mathematically  by  notinq 
our  sample  set  wj  should  allow  us  to  estimate  the  signal  which  created  the-  or 
(drop  ?-  weighting) 


f(t) 


F(,Ve 


-ico.t 

J 


(1 


We  have  reconstructed  the  source  of  the  original  data  and  can  use  it  in 

tins')  to  f  i  n  ri 
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In  principle,  the  estimate  of  F(  is  a  weighted  sum  cf  all  F(w.)  samples  and  on 
the  order  of  N  operations  arp  indicated  for  each  desired  data  point.  As  a  prar.tica1 
matter  we  see  that  nnlv  those  samples  within  a  frequency  space  resolution  piemen* 
contribute  significantly  to  the  estimate.  For  adequate  samplinq  we  would  have  two 
such  points.  It  may  he  nice  to  add  more  by  oversampling  at  the  video  and  let.  ‘he 
i nterpol a* ' on  serve  a  pres  urn,  as  well  as  a  formatter  function.  Practical 
implementations  obtain  good  estimates  with  only  a  near-neighborhood  data  set.  Fer 
two-dimensional  data,  the  indicated  interpolation  process  is 


■V  V  1  oinc  (i"V  l sjnc  (6"V  2 


with  N"  operations  of  each  of  N  points  in  our  minimally  sampled  case. 

In  one  dimension  we  may  alternately  choose  to  implement  an  interpolation  process 
based  on  a  Taylor  series  representati ve  of  reasonably  behaved  data.  We  have  date  a  + 
position  * (<  -  nh  whore  n  is  an  integer.  We  desire  f(x  *  A)  where  we  can 
reasonably  assume  .  I'sing  a  Taylor  series  we  find 


’  x  )  +  f 1 


mi  roqij-irn  estimates  n*  fl,f",f"',  efr,  at  x.  Treating  the  available  data  in  a 
TAv’’r>r  or  r-ifj-jtp  difference  manner  we  have 


fly  •  nh)  f'xi  *  nh  f1  *  4  • L  f'\  + 

x  2 '  1  x 


It  is  possible  to  write  a  series  of  finite-difference  equations  usina  MP8'  arvJ 
obtain  estimates  of  the  derivatives  of  f  at  x  with  pven  more  accuracy  as  more  o*  *he 
neighborhood  samples  are  involved.  If  a  first-  ,md  second- order  ost  irate 


a  de  n -ate.  firm 


•  ■*  ■'  - 


"-“.vZvTy v  ‘ 


f  s 
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« 


f(x+h)  -  f  (x-h) 
2h 


f(x+h)  +  f(x-h)  -2f(x) 
.  2 


We  can  proceed  to  estimates  of  higher  order  derivatives  by  considering  more  equations 
and  sample  points.  The  extension  to  two  dimensions  falls  apart  for  the  radar  case  as 
regular  samples  are  generally  not  available  in  the  orthogonal  direction.  We  could 
associate  a  unique  h  with  each  data  point,  but  the  forms  become  more  complex.  The 
need  for  2-D  interpolations,  i.e.,  consideration  of  diagonal  elements  vs  sequential 
1-0  interpolations  to  estimate  a  data  point  value,  is  one  of  accuracy.  The  selectior 
seems  a  matter  of  choice,  and  on  the  surface,  appears  to  be  a  tradeoff  between  a  near 
neighborhood  and  somewhat  larger  linear  neighborhoods. 


The  above  is  an  attempt  to  illustrate  that  a  data  set  belongino  to  a  wel1  behaved 
process  can  be  interpolated  to  another  data  set  representing  the  same  process.  How 
this  is  done  is  not  of  particular  import  in  this  overview  of  the  concepts  of  SAR. 

All  we  need  to  note  is  that  our  ideally  formatted  data  indicated  a  minimal 
requirement  for 


N  1  og^ 


BBSS 


■-  .VJV  _>  A  -  »  .  ■  .  -  .  .  A  .  -  .v  A  .  .  '■ 


<■,  A  A  ■ . 

■  -  ■  •  r  ' ' y. 


‘  ipf,  an  .  v  N  71  x  e  1  s  inaqe.  We  il  Uist  rated  above  narh  r 

P-u  i  ,-ed  V  ’  r*pl;f-  estimates  1  irynl  ve  some  number  nf  opera*  ’  or-  s  .  To 

pf :  ,  Mr,.'/  *  >'••>  •  "t^rpr. 1  ^  np  loan,  wp  would  have  to  examine  a  »»' ver  pr*-- 

n  .'Ip*.)’1.  f  >  ■  w ,  wp  r  '  hote  *hat  fhp  interpolation  process  dominates 

o 

*.  r .  "'oo '•  •  •-;•• 1  4  *>.  for  each  of  the  M  -  data  points  exceeds  the  cor:p1ov 

■■q,\  Hi.’*  t  r-  ’ps.  ’’f't  .  ;>’e>  *.  r;.  v  bp  some  conbinations  0*  hardwa  rr  arm  sot 

rtp>'pp'' c>*  inn  n.v,  .'f  r  *■  be  *rp  dominating  process  in  some  recen*  SAP 

r:v  nnen*;.  *  :  a  r  •  a  mot  ••H  •  *  ud;es  fo  reduce  its  intensity  ,  or  ideally  tt 

1  im-' 0  a  *  e  *  tv  at  i  v»  j-r,  mssino  schemes  which  will  not  s  ianifican* 

2 

rr-c'-e  the  '  1  1  :  *=-,•.  -.iCn  op./  -nd  *he  equivalent  of  N  log„M  Putter*]  -e 


r.wa  r^ 
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PACK-PRO, iFCTION  PROCFSSING 

This  section  is  included  to  illustrate  a  variation  which  mav  ho  ^nplemen*o(<  " 
process  the  data.  The  appendix  details  another.  They  are  offered  as  i  1 1  uc  *  ?-<i  *  >  ops 
and  no+  necessarily  proposed  as  having  merit.  Rack-projection  prnressinq  i r- 
basically  a  piecewise  evaluation  of  an  integral  or  summation.  For  rhp  cases  of 
interest-,  we  can  consider  the  1-P  or  2-P  discrete  inverse  Fourier  transforms.  In 
1-P  case  we  can  break  the  summation  of  (93)  into  a  series  of  subintervals  and 
estimate  the  desired  quantity  as 

S  (kAu)  =  E  S  (kAuj)  'i  - 

o  o  m  1  - 

m 

wherp  \  s  :n  <  N.  The  quantities  summed  are  the  results  of  surnmino  suMntervals  if 
the  range  of  j  or 

j+1-1  •  9 

n'  -l  2tt  . 

SQ(kAuj)m  =  E  S.(kAw)e  N  11  • 

where  we  designate  the  onset  of  a  subinterval  by  j  .  It  should  be  obvious  that  we 
can  do  this;  however,  the  implementation  is  different.  Implementing  (93l  usino  ar 
FFT  algorithm  results  in  a  rippling  data-reduction  sequence  with  the  estimates 
(k  -  w)  appearing  out  of  separate  spigots  some  time  later.  The  form  of  (1131 
illustrates  that  each  subinterval  provides  a  complex  increment  for  our  estimator  o* 

■>o  (  k'.w),  as  did  each  term  of  (93).  For  each  of  the  k  aw  of  interest  a  data  store 
has  to  be  established.  The  results  of  each  subinterval  summation  is  then  added  to 
the  quantity  already  in  store  as  a  result  of  past  subintervals  summations.  When  n’i 
samples  have  been  treated  once,  and  only  once,  this  sequential  summation  will 

produce  the  desired  estimate  (112).  There  are  no  restrictions  on  sample  lenqths. 

The  extension  t.o  two  dimensions  proceeds  in  the  same  manner;  only  now,  +  he 
subintervals  can  be  viewed  as  subareas  of  the  total  area  of  the  data  space.  In  this 
'•ase  we  reouire  a  data  store  for  each  pixel  in  the  final  imaoe  array,  q  *  •,  for  what 

we  have  discussed  above.  The  data  aperture  can  be  broken  info  as  many  nonoverl appi rn 

subareas  as  desired.  The  item  of  note  is  that  the  data  store  and  its  implied  '  •  . 


T  P  _  Q  h  - \ Ifp 


''pipy  'Mjphf  f  *■  O  one  .('rp  I,1 


'  n  ♦  he  PAD  irpns  ,  a  nne-d  i  mer s  i  or>  ■»  ’  p» ; r  a  1  spatial  1  ’  ph  *  •’fidi,'  •  ♦ 

•♦spI*  as  a  h  vh  r*  •»  d  hack-pro ieotor  svstfv'.  ;.e,,  tcnus  to-npt  i  r  Price  *■’'  ’  >  , 

■■  ,'1-,r'.-pp*’ija  1  ■ v  ^r-asihlp  to  reduce  the  :  '  .r-t  ormaf  ted-SAP  da 4  a  '•  !  >  -  . 

:  :~r.c  c  ■>  ncs  nf  thp  data  i  n  s'!  ’fjhtl  v  we-"1  eg  su  hare  as  c*  ♦ho  da  t  a  s>  v  e  «  ;• 

*  ’  '  '.s-ratpd  in  Diqurp  11.  !n  practice,  those  si  ices  would  hP  arprnpr ' a  * n 1 ,  wt  ■  * 

sa"’p’es  n  +"  «-hp  da  fa  for  any  --  .  rnr  a  oh  i  rp  inpleniPntat  ioti ,  the  da  ♦  a  w  u 1  d  h"  ’  •  • 
Kserha^d  v’dpn.  o*  anv  radar  PRF  line  superimposed  on  a  carrier  to  rnpj>--*n  *t< 

’i.j':a*e  optical  imaoes  produced  by  the  Bra  no  cell  wher  i  1 1  unina  *r-d  tv  a  rchc-.  >  * 
v.o-,f  *ront.  the  complication  tor  the  optical  processing  is  *hat  a  coherent  *  .• 

v'mv  r-f  t,  x  ^  e  1  erne  n  t  s  has  to  be  provided.  A  spatially  and  tpmporar ’  ’  v  'ohor,-..  • 
-nS-rpn'P  wavetonnt  a^o  has  to  be  introduced  before  detection  *o  pv*<-.r+  h-*h 
•~r’i tgde  and  phase  nt  the  signal  on  each  detector,  tor  each  subapnr*|jre.  Thir 
task  .  "hose  amplitude  and  phase  data  are  then  supplied  *n  the  ti.rrinq  s*o,-, 

•  -  - ’  ■> ted  w^h  each  detector. 

Tho  issues  w’th  hack-prniect i on  proenssinqs  should  nn+  he  "are  +hev  are 
pns s i h 1 e ° "  the  iesups  are  engineering  attributes  ot  anv  implementation.  A  ma > 

'  s  s 1 1 r  ’  s  the  noice  introduced  bv  errors  in  the  process;  in  particular,  noise  cr 

C 

eytr,iepous  siqnals  which  tend  to  be  coherent  and  limit  the  potential  ’/„  or  dru 
’••Jr.;es  «t  the  proressed  data.  An  example  ot  such  a  coherent  noise  mav  he  a  rippi.. 
a  power  supply  common  to  many  elements  nf  a  phased  array,  when  components  ot  earh 
: r r a y  element  have  a  pa  +  h  delay  which  is  voltage  sensitive.  The  oues+  ions  whir1"  t. 
*c  he  addressed  in  back-proiection  schemes  center  around  how  well  a  practical 
■rp'onppt  approaches  +he  ideal  ot  the  iterative  summation  nt  QlZ)  or  its  2-H 
counterpart.  These  topics  have  their  counterpart,  in  digital  processings  where 
quart ’ zat i on  levels  and  tunings  have  to  he  addressed. 


^ ^ I**—.  W"^. al" «.  m  W V-  f  _ V-  ■/_ 
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'VsCRSSION 


The  preceding  discussion  has  separated  t.he  basic  physical  requ  i  remen  t  s  •*  C'A 
’run-  the  engineering  concerns  which  arise  during  any  implementation  of  tre  ire 
We  ^ee1  that  the  physics  of  such  a  radar  are  those  of  a  microwave  inter*  emmo* 
matched  filter  which  can  probe  the  data  generated  for  the  presence  of  a  <rv— ■ 
sianal.  The  implementation  requirements  of  such  an  interferometer  are  t-as1', 
specified  and  are  built  into  any  SAR.  There  are  no  basic  restrictions  on  *  •-,>>-cr 

waveform  or  platform  motion.  The  variations  and  activity  in  SAR  renters  n;  >r 

implementation  of  the  matched  filter,  and  it  has  at  least  been  implied  that  t h(.< 

efforts  are  basic  to  a  SAR.  They  are  not. 


*  t 


V’ewed  with  hind-sight,  the  development  of  SAR  processing  techniques  seems  to 
have  evolved  more  along  the  lines  of  available  techniques  and  the  means  to  over'"'"-' 
‘he  limitations  of  previous  techniques.  Initially  all  SAR  data  was  optica1’'' 
processed  and  a  range-azimuth  data  format  was  perhaps  rather  obvious.  As  dMv'ta1 
processina  was  applied  to  strip-mapping  SAR 1 s  the  range-azimuth  format  remav-ed. 
however,  the  azimuth  data  in  any  range  bin  was  prefiltered  so  a  batch  process  rop'd 
he  used  to  simultaneously  extract  the  data  on  all  azimuth  cells  in  a  given  range  b->. 
'his  afforded  a  computational  efficiency  since  the  same  data  did  not  have  ‘o  kp 
in  many  processings  as  it  is  used  in  the  optical  processor,  or  as  it  would  be  usei  ;> 

:  brute-force  matched-filter  implementation. 

par,ge-azimuth  formatting  has  remained  the  mainstay  of  SAR  processing  roncep+s. 

Put  imposes  rather  severe  limitations  on  the  data  space  which  can  he  processed  t  • 
ha*<-h  mode  as  resolutions  improve.  These  limitations  were  ^irst  encountered  d  i  j  '  <  -  ■  : 
experiments  with  the  rotary  platform  SAR  emulator  with  its  fine  resolution 
capability.  Attempts  to  overcome  these  limitations  prompted  the  polar  format 
■'hserva t ions  of  Walker.  In  the  polar  formatting  procedure  the  matched  filter  f»r  - 
reference  position  is  removed  from  all  data  and  the  residue  is  concept ua 1 1  y  1 

• r  a  po 1  a r  forma1'  for  additional  batch  processing.  During  the  rejection  ct  on ‘a  1 
this  ffnrat,  individual  scatter  points  may  move  throuqh  many  resolution  and  i 

range-azimuth  format  becomes  impossible.  For  points  in  the  near  neighborhood  of  *hn 
r(  'nrnnrp  position,  the  frequency  aspect-angle  polar-format  ted  residue  h^s  tgP 
f  sea  * ‘er-pns i t i on-un ioue  ?- D  spatial  monotones.  The  process i  ru  n‘  tho?r  - onot  nr.  • 

’’ 1  ‘  -1  spectral  analysis  of  the  data  which  exists  in  some  ,ip<T*  u*'<'  r*  *  ho  f 1  : 
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space.  Much  of  the  present  activity  in  Spotlight  SAR  centers  on  the  digital 
processing  of  this  polar  formatted  data.  The  spectral  analysis  is  rather  nicely 
accomplished  using  implementations  of  the  FFT  algorithm.  However,  this  alnorithm 
comes  with  the  burden  of  a  particular  sample  length  and  regular  data  samples.  The 
FF^'s  need  for  regular  data  samples  forces  an  interpolation  from  measured  valuer 
along  polar  lines  to  positions  appropriate  for  the  FFT  algorithm  processes.  Th->'s 
interpolation  tends  to  be  more  computationa 1 ly  intensive  than  the  final  data 
^eductions  with  FFT's  and  has  been  the  focus  of  considerable  activity. 

While  the  present  focus  on  polar  formatting,  interpolation,  and  FFT  algorithms  is 
very  well  based,  the  purpose  of  this  report  is  to  illustrate  that  these  are  not 
fundamental  radar  processes.  They  are  processes  associated  with  a  particular 
approximation  to  the  ideal  and  as  such  come  with  their  new  set  of  limitations.  These 
1  imitations  are  treated  elsewhere.  We  feel  that  the  report  illustrates  this  position 
ang  a’so  describes  the  basic  characteristics  involved  in  various  implementation 
forms . 

With  the  background  and  illustrations  provided,  a  reader  should  be  able  to 
separate  basic  physical  requirements  from  engineering  considerations.  The  concepts 
n*  polar  formatting  are  very  powerful,  and  the  use  of  FFT's  to  reduce  the  residual 
data  may  well  be  optimal  in  many  situations  where  FFT's  and  their  overhead  are  not 
optimal.  The  concepts,  approaches,  and  procedures  of  this  report  may  help  to 
objectively  access  how  best  to  approach  the  processing  related  aspects  of  a  desired 
SAR  implementation  and  not  have  to  blindly  accept  approaches  taken  in  past 
implementations,  assuming  them  to  be  inviolate. 
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APPENDIX 

Segmentation  of  the  Discrete  Fourier  Transform 

In  this  Appendix  we  will  look  at  a  segmentation  approach  to  the  discrete 
transform  illustrated  by  (93).  As  in  (93)  an  assumption  of  regular  sampling  is 
invoked.  We  will  address  the  digital  analogs  of  the  multiple  local  oscillator 
radar,  discussed  in  the  section  on  Processing  Considerations,  which  come  abou1  ’;r  a 
attempt  to  introduce  batch  processing  to  the  range  data. 

Let  us  go  by  (93)  and  consider  that  for  some  reason  we  are  interested  in  only  f’ 
of  the  N  possible  estimates  of  S  (k  Aw)  given 

N-l  -  i  2_tt 

S  (kAw)  =  I  S.(kAx)e  N  J  s.i 

j=0  1 

the  above  represents  an  optimal  sampling  of  the  data  space  centered  on  zero  freguen 
and  an  optimal  sampling  of  frequency  space  to  extract  all  significant  data.  Let  us 
begin  by  assuming  M  is  a  power  of  2  and  the  band  of  frequencies  of  interest  are 
centered  around  v  =  0.  The  process  indicates  M  x  N  brute-force  operations  to 
accomplish  the  desired  estimates. 

How  do  we  reduce  this  since  we  have  an  excess  of  data  points  for  the  band  of 
'^equencies  of  interest,  but  our  sampling  is  driven  by  the  input  data?  Tn  our 
previous  discussions,  fhe  analog  sections  of  the  radar  provided  IF  filters  and 
eliminated  the  uninteresting  signals  before  we  considered  the  frequencies  of 
interest.  Filtering  can  also  be  implemented  digitally,  and  we  can  consider  the 
rpsult  of  summing  b  adjacent  samples  from  a  monotonic  data  set  with  a  basic  signal 
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B  .  . 

2  hertz 


Sampling  at  a  rate  a  =  and  considerinq  h  samples  beginninq  at  ti>p  *=  * 

t  B 

arrive  at  a  composite  sum 


S 


bv]  e  i4t0tmAt] 
0 


Performing  the  indicated  summation  yields 


b  i  sin  ("?->  « 

e  ‘  T1  ^  - - \7~~~ 

Sin  2  At 


Note  that  S'  appears  as  a  weighted  sample  of  the  original  signal  at  time 
t  =  tQ  +  (b  -  1 )  a  ^ .  The  weighting 

sin  (^)At 
sin  (^)At 

is  the  array  factor  seen  wherever  multiple  coherence  siqnals  interface,  ’.e., 
collections  of  microwave  or  optical  sources  which  are  arrayed  in  space  and/or 
frequency.  The  weighting  is  frequency  dependent  with  peak  amplitudes  of  b  at 


For  T  -  1/B  we  find  peaks  at  mB.  Our  only  interest  is  with  the  peak  at  w 
as  our  maximum  frequency  is  +  R/2.  The  3-db  width  of  the  peak  is  approximate! y 
but  the  function  (A-6)  falls  off  relatively  slowly  and  has  significant  sidelobf>c 
beyond  B/b  which  may  impact  downstream  processes. 
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To  overcome  the  filter  sharpness  and  sidelobe  influence  the  samples  can  he 
weighted  in  amplitude  before  summation  (Taylor,  cosine  squared,  binomial, 
Tchebvscheff ,  etc.).  The  result  of  these  weightings  is  a  sharpening  o*  the  r 
of  the  array  factor  (analog  of  (A-6)  for  the  particular  weighting  i  and  a  virtu*' 
elimination  of  the  sidelobes.  These  advantages  come  at  the  expense  of  a  suppre'a- e'o 
of  the  peak  amplitude  (w  =  0)  by  perhaps  l-?db  and  a  sliohf  broadening  nf  +he 
bandpass  to  something  slightly  greater  than  B/b.  He  have  created  a  filter  which 
preserves  the  strength  of  the  signals  in  some  band  approximatina  B/b  and  have 
destructively  interfered  those  outside  this  band. 

For  our  case,  (Al),  we  may  want  to  do  somewhat  more  than  1  of  these  presums 
alnna  the  data  record  length  to  provide  data  for  our  estimation  of  Sjl  '  w).  Thi 
would  be  oversampling  so  that  our  estimates  of  the  samples  of  interest  would  no* 
have  members  too  close  to  the  skirts  of  the  filter.  To  accomplish  this,  the  blocks 

of  data  which  are  presumed,  are  separated  by  less  than  b  A*  ,  i.e.,  leading  and 

t 

••railing  samples  are  components  of  two  presums.  The  resultant  M  -1  samples  may  row 
bp  processed  using  an  M  +  point  FFT  which  requires  on  the  order  of 

M  log2  M 

operations.  This  may  appear  significantly  less  than  the  MN  brute-force  operations 
indicated.  However,  we  have  to  include  the  presuming  operations  to  see  if  wo  have 
gained  anything.  If  we  consider  weighting  and/or  frequency  shift,  to  baseband  (in 
oeneral  our  subband  may  be  anywhere  in  the  total  bandwidth)  we  will  require  on  the 
o*der  of  Mb  preconditioning  or  filtering  operations  for  a  total  process  approximate 


M  [log2  M  +  b] 

operations.  As  *1  Mb  we  may  write  (A-9)  in  the  approximate  form 

log.,N  -  log,,  b 

b 

The  above  indicates  a  potential  for  significant  savings  over  a  brute-force  estimate 
of  the  M  values  of  S  (k.tw'  which  are  of  interest.  However,  if  our  intprer*  "c 
really  with  the  full  N  values  and  the  filtering  was  only  an  ut’erpt  to  ■ 
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analog  function*;  of  the  radar,  we  lose.  We  would  require  b  such  pnncpdures  for  +  hn 
N  total  points  of  ktw.  From  A-10  we  see  the  totality  of  operations  would  export 
*tp  N  loq^N  FFT  operation*  indicated  for  the  original  data  set. 

This  appendix  has  been  included  to  provide  a  brief  introduction  to  digital 
‘ilferino  and  to  further  illustrate  the  power  of  efficient  data  reduction  a i qori f hms . 
F'n  the  surface  it  appears  we  would  like  to  condition  the  input  da  t  =  tor  the  larqest 
possible  FFT,  or  the  total  span  of  interest  if  wp  can  build  an  arbi trary-sized  FFT. 

Tf  only  a  small  portion  of  the  collected  data  is  of  interest,  t-ntprinq  can  provide 
efficiencies.  A  tradeoff  may  be  possible  if  considerable  conditionina  has  tn  be 
done  on  the  raw  signal,  and  the  filtering  operations  can  be  incorporated  into  this 
process  while  incurring  little  addition  complexity.  We  could  be  left  with  an  N 
point  FFT  with  its  N  log^N  Butterflies,  or  b,  M  point  FFT's  with  N  log?M 
Butterflies  where  N  =  Mb.  This  is  an  analog  of  the  radar  case  where  we  ignored  the 
increased  complexity  of  the  analog  sections  of  the  radar  and  found  our  processing  wa* 
mirimized  using  the  shortest  possible  FFT. 
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